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Abstract 


Laser-based spectrophotometric methods which have been proposed for 
the detection of trace concentrations of gaseous contaminants include 
Raman backscatteririg (LIDAR) and passive radiometry (LOPAIR) . In this 
report, we discuss remote sensing techniques using=.lase'r^spectrometry , 
and in particular a simple long-path laser absorption method (LOLA),, 
which is capable of resolving complex mixtures of closely related trace 
contaminants at ppm levels. 

A number of species were selected for study which are representative 
of those most likely to accumulate in closed environments , such as submarines 
or long-duration manned space flights. Absorption coefficients at CO 2 laser 
wavelengths were measured to an accuracy of + 1% or better, for each of 
these species. This data base was then used to determine the presence and 
concentration of the contaminants in prepared mixtures of twelve to fifteen 
gases. Computer programs have been developed which will permit a real-time 
analysis of the monitored atmosphere. Minimum detectable concentrations for 

Individual species are generally in the ppm range, and are not seriously 

\ 

degraded by interferences even in complex mixtures. Estimates of the dynamic 
range of this monitoring technique for various system configurations, and 
comparison with other methods of analysis, are also given. 
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CURRENT MONITORING TECHNIQUES 


A. Backgrotmd 

Pollution is not merely the addition of substances to the air as a 
result of man's activities. Most of the substances considered to be 
pollutants are already present In the environment from many natural 
processes. There does exist, however, an assimilation capacity of the 
envlronn^nt (ACE) , ^ which through dilution Is able to maintain a 
balance. The ACE Is fairly substantial for most substances but not 
Infinite. If the addition of substances to the environment Is below 
the ACE, the effect Is not noticeable, and no "pollution" results. In 
the last century the growth In population, along with a corresponding 
Increase In wealth and consumption, has resulted In Increased waste per 
capita, which Is severely straining the ACE, especially In urban areas. 
Common property resources (air, water, public recreational facilities) 
are saturated. In addition, present Industrial technology Is moving 
toward use of non-blodegradable substances (ACE Is very low or zero for 
aluminum beer cans) or substances like pesticides which are cxnnulatlve 
In organisms high In food chains. 

The air and water are viewed as free, and consequently these resources 

are wastefully used. Until some value Is placed on these environmental 

2 

resources, the present misuse will continue. The social costs of 

• > ■ ' ’ 

pollution Include maintenance and cleaning costs, loss of recreational 

facilities, and the ultimate consideration, effects on health and 
mortality. These will be discussed further below. 

Changes can only be brought about by placing a value on clean air. 
Suggested methods Include: 
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1. Incentive pricing by means of graduated fines, redistribu- 
tion of the burden of paying for pollution, and adjustments to receptors; 

2. Federal subsidizing to help industry change generating proc-. 
esses, encourage recycling, and to provide tax incentives for complying 
corporations; and 

3. Regulation through use of large flnes^’luid strict non-beneflclal 
enforcement . 

Regional management is needed for complete control of problems. A program 
combining all three approaches would be best. An Integral part of such 
a program would be an accurate monitoring scheme both to force compliance 
to the law and to permit the study of pollution dispersion and thus the 
program's effectiveness. 

3— ] 

The sources of various atmospheric pollutants have been listed often. 
However, their physiological effects have not received as much publicity. 

In order to form a basis for understanding the toxicity of polluted atmos- 
pheres, we will consider the effect of pollution on human health next. 

Effects of individual pollutants are relatively easy to assign. Carbon 
monoxide bonds with hemoglobin 210 times more readily than 0^, and thus 
requires lungs and heart to work faster.^ At levels presently found in 
downtown areas, sensory response is io^alred. At higher concentrations 
symptoms include headache and nausea. There appears to be no threshold 
below which symptoms are not in evidence. SO 2 affects the upper 
respiratory tract, constricting airways and slowing mucous flow, permitting 
bacteria to have longer residence time in the lung. NO^ destroys alveoli 
(which permit gaseous exchange between the blood and air) and destroys 
lower lung cilia, which sweep bacteria and particles out of the lungs, 
thus lowering resistance to pneumonia. A secondary effect of NO 2 (and 
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chlorofluoromethanes) pollution night be the destruction of the ozone 

layer, ^ permitting solar ultraviolet radiation . to reach the surface of 

the earth with a resultant Increase In the incidence of skin cancer. 

Chlorofliioromethanes are relatively non-toxic and are widely used as 

refrigerants and propellants. It must be said, however, that the 

actual effect of chlorofluoromethanes on the upper atmosphere Is still 

a matter of considerable uncertainty and controversy. Ozone disrupts 

the cell membranes, causing fluid leaks Into the lungs, and lowers 

the diffusion of gases In the alveoli. At high concentrations, headaches 

and lethargy are symptoms. Fortunately, 0^ concentrations Inside btjlldlngs 

are usiially much lower than outdoors, due to the rapid first-order decay 

8 

of 0^ on rubber, plastic, or textile surfaces. 

Hydrocarbons generally are not toxic In the concentrations present 
In most polluted atmospheres, but participate as reactants In photo- 

9 

chemical smog, forming peroxyacylnltrates (PAN) and other carclnogeiu. 

Some widely used commercial substances are also suspected of causing 

cancer. These include CCl^ (liver, kidney), CHCl^ (heart, liver, kidney), 

and vinyl chloride (liver, kidney, lung, spleen, bone marrow) . 

The threshold limiting values (TLV's) set for several common pollutants 

for work room air are listed In Table I. Table II shows concentrations 

12 

found In "clean" air In urban environments and the level set by the 

13 

Environmental Protection Agency (EPA) In 1971. Many feel that the EPA 
permitted levels are too high. A consideration of the levels during 
past pollution episodes will provide a perspective. 

In London in 1952, a four-day fog claimed an Increased fatality 
count of over 4,000 people, mostly the old and the very young. PoUuclon. 
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Table I 

Some Maximum Permissible Contaminant Levels In Closed Environments 


Compound Work Room Air Ambient Level Extended. 

(8 hour exposure In ppm^® Period - 6 months^O in ppm 


acetaldehyde 

100 

10 

acetone 

1000 

300 

acetonltrite 

40 

1 

ammonia 

25 

25 

benzene 

10 

1 

t-butanol 

100 

40 

2-butenone 

200 

20 

carbon monoxide 

50 

15 

cyclohexane 

300 

60 

o- dl chlorob enzene 

50 

5 

Freon 12 

1000 

100 

Freon 113 


50 

1 , 2-dichloroe thane 

50 

10 

ethyl acetate 

400 

50 

ethanol 

1000 

50 

furan 


0.04 

isopropanol" 

400 

40 

methanol 

200 

40 

methylchloroform 

350 

50 

nitric oxide 

25 


nltrogendloxlde 

5 

0.5 

ozone 

0.1 


sulfur dioxide 

5 

1 
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Table II 

Pollutant Concentrations In Atmosphere In ppm 


Conq)ound 

Background^^ 

12 

Orb an 

(annual mean) 

13 

EPA Standards* 

SO 2 

2-4 X 10“^ 

3 X 10“2 

3 X 

10~2 annual mean 




1.4 

X 10-1 24 hour max. 

NO 2 

1-3 X 10“3 

5 X 10"^ 

5 X 

10“2 

CO 

1 X 10-1 

4 

9 

8 hour max. 




35 

1 hour max. 

O3 ^ 

1 X 10“2 

3 X 10-2 

8 X 

10"2 1 hour max- 

CH4 

1-1.5 

2 

— 


Hydrocarbons 
other than CH^ 

1 X 10-3 

5 X 10"1 

2.4 

X 10“^ 3 hour max. 

* The 

maximum values are 

permitted only 

once 

each year. 


in London Is primarily SO 2 based (stationary source combustion)^ and 
yet the highest dally averages were 1.35 ppm (parts per million) SO 2 
on the last two days of the fog. This level Is only ten times the EPA 
maximum level, and thousands of deaths occurred. In New York City In 
1966, a three-day eposlde resulted In 168 extra deaths and the hourly 
SO 2 levels remained between 0.3 and 0.7 ppm with only two readings at 
1.0 ppm. This episode had pollution In excess of EPA standards by a 
factor of only 3 or 4. Although some of the people who died during 
these acute episodes may haye been more sensitive than the general 
population, "It seems more likely that they succunbed to the combined 
effects of sulfur dioxide and other pollutants that simultaneously 
accumulated during the episode. This Is further substantiated by 
the Donora, Pennsylvania Incident In 1948 when 40% of the populace 
was affected and 10% had severe respiratory difficulties that per- 
sisted for up to a year. Some of the fatalities had no history of 
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chronic disease. It Is fair to state that smokers and those with a 
chronic cough, as well as the very young, are likely to suffer most 
during these episodes. 

Controlled studies on the effects of mixtures of air pollutants are 
rare, but In general show a reduced tolerance, especially among the 
sensitive. Consider an example of ctmmlatlve effects: NO2 Irritates 

the lungs, reducing pulmonary efficiency, O3 Interferes with the trans- 
port of O2 from the blood to the tissues, SO2 Increases airways' 
resistance, hindering breaching, CO binds with the blood preventing 
O2 transfer, and Pb prevents the synthesis of new blood, aill resulting 
In a serious oxygen deprivation problem. Thus, realistic standards for 
mixtures cannot be set by the current technique of partial threshold,^® 
and thorough Investigation of possible synergistic effects of mixtures 
should be carried out. However poorly defined these standards, pollution 
must be held to these levels. The resulting need for monitoring devices 
for atmospheric study, for enforcement, and as a warning device provides an 
Inqietus for much current research. 

B. Monitoring Networks 

Although many types of monitoring schemes have been proposed, they 
may all be classified as either point sampling monitors or remote monitors. 
To date aliujst all data gathered about ambient air quality has been through 
the use of point sampling networks. This type of monitoring network pro- 
vides information only about the conditions In the Immediate vicinity of 
the sanq>ler. Strategies must then be devised to gxiess conditions at all 
other locations between sampling stations. Obviously, point sampling 
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data may provide an accurate measure of the ambient pollution In a given 
area, especially In submarines and future spacecraft where good central 
ventilation performs thorough mixing and to a lesser extent In factories 
and tunnels, where there Is good Interior mixing, but poor exchange for 
fresh air. Point sampling of specific sources, such as smokestacks. Is 
another Important application. The capability for remote monitoring 
adds an additional dimension to a network. It can provide air quality 
Information at many specific urban areas or sources, as well as vertical 
pollution profiles and averaged pollution concentration over general 
downtown areas. This Information Is necessary to validate pollution 
dispersion models and the effects of Individual sources on general 
air quality. The remote techniques all rely on the Interaction of 
radiation with the pollutants. 

Point sampling techniques lised situ provide Information of the 
air quality where people are working or spending a lot of time. It is 
relatively Inexpensive to collect remote samples (In canisters, bottles, 
or bags) and analyze them all at a central location, using accurate 
sensitive techniques capable of detecting a variety of contaminants. 

Such a program would be unable to provide real-time Information about 
air quality, thus providing no warning when emergency situations exist. 

A remote monitoring network would not suffer this limitation because- the 
signal may be processed as soon as it Is received. Unfortunately, unlike 
the point sampling schemes for which widely applicable analyzers exist, 
most remote methods at present have been able to provide Information 
about only a few species, and sensitivities and reliability have not yet 
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reached the level attained by the best point sampling techniques. However 
these limitations are technological » not fundamental. In nature and the 
present experimental work Is directed toward solution of these short- 
comings. 

Real-time Information Is vitally liiq>ortant in the closed environments 
found on long-duration manned missions, as In spacecraft or submarines, 
where sudden contaminant level changes due to leaks or breaks are 
quickly spread through the entire environment. In addition, because the 
concentration of substances which are ordinarily present In harmless trace 
amounts can slowly accumulate to toxic levels, the monitoring system must 
be able to simultaneously detect levels of many contaminants over a 
range of concentrations. 

In these closed systems where the primary sources of contaminants 
are outgasslng of non-metals, metabolic b 3 rproducts, and fluid leaks, 
the problem of altered toxicity due to the presence of many other contami- 
nants Is serious. Present point sampling monitors for these environments 
Involve the use of charcoal canisters with analysis following the mission. 

Canisters flown aboard the Apollo and Skylab missions yielded over 300 

16-18 

co&q>ovmds upon desorption. The quantities of some typical sub- 

stances fo\md In the canisters are given In Table III. However, there 
exists no simple relation between these measurements and the actual 
ambient level of the compoimd at any given time; even the relative values 
are misleading due to differing adsorption coefficients' for the different 
compotmds. A striking example of the need for real-time Information was 

the contaminant leak on descent during the recent Apollo-Soyuz mission, 

19 

where Ignorance of the irritant hindered medical treatment. 
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Table III 


Maximum Contandnant Level Found upon Analysis In any 
Canister on Manned Missions, In yg/g Charcoal 


Contaminant 

VII 

Apollo 
VIII IX 

Missions 

X 

XI 

XII 

Skylab Missions 
3 mid 3 end 4 mid 

acetonitrile 

0.97 

— 

0.15 

0.33 

0.78 

■ 

0.0 

0.0 

0.20 

benzene 

0.6 

3 

6.19 

2.29 

5.59 

5.06 

0.0018 

0.072 

0.0068 

tert-butanol 

0.61 

2.2 

0.18 

0.00 

0.079 

1.2 

0.0 

0.0 

0.0 

cyclohexane 

14.3 

0.5 

5.35 

2.27 

4.82 

5.67 

0.11 

0.075 

0.03 

1,2 dichloro- 
ethane 

— 

0.0 

— 

— 

— 

— 

0.0 

0.0 

0.0 

o-di chloro- 
benzene 

0.008 

0.02 

0.002 

0.0012 

— 

— 

0.0 

0.0 

0.0 

ethyl acetate 

8.3 

17 

6.22 

2.47 

4.87 

4.56 

5.5 

3.4 

0.20 

Freon 12 

125* 

83* 

412* 

91.3* 

3.03 

1.02 

0.30 

0.49 

0.12 

Freon 113 

1300 

328 

276 

353 

160 

214 

3.6 

0.97 

0.28 

furan 

— 

0.8 

0.57 

0.33 

1.04 

1.31 

1.3 

0.075 

0.058 

Isopropanol 

19.2 

3.8 

30.2 

3.4 

1.06 

19.2 

61. 

17. 

12'. 

methyl 

chloroform 

4.55 

2.0 

0.66 

0.31 

5.89 

602. 

— 

— 

. 

methylethyl- 

ketone 

7.2 

5.4 

1.7 

0.16 

0.25 

1.74 

0.0 

0.0 

0.0 

vinyl 

chloride 

0.13 

— 

; 

— 

— 

— 

0.0 

0.0 

0.0 


* Figures for Freons 12 and 22. 

Freon 12 is dichlorodlfluoromethane; Freon 113 Is trichlorotrlfluoro- 
ethane; Freon 22 is chlorodlfluoromethane. 
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C. Point Monitoring Techniques 

The original point sampling standard techniques Involved stoichiometric 

reaction of preconcentrated samples of each pollutant, and although these 

methods are very accurate and sensitive, they require many reagents, a 

trained technician for analysis, and a long collection period of one 

hour to one day, thus destroying real-time fluctuations. There has 

been a constant trend toward more mechanized analysis methods requiring 

12 


fewer support chemicals and allowing more routine analysis. 

Gas chromatography offers both sensitivity and the ability to monitor 

a great many species after calibration. Detection limits are In the ppm 

to sub-ppb (parts per billion) range, depending on the detector iised at 

4 21 

the output to the column. * Although a carrier gas and preconcentration 
are needed, analysis Is routine for expected compounds, and this technique 
is vastly superior to the wet chemical methods previously used. Urban 
real-time observations are further prevented by the necessity of collecting 
samples peripherally for -later central analysis. It would not be economi- 
cally feasible to locate a gas chromatograph at every monitoring station. 
Presently, non-unlformltles, degradation of the column material, and 
clogging of the preconcentration beds, as well as the necessity of complex 
flow controls, diminish this system's reliability and reduce Its desira- 
bility for submarines and spacecraft. The need for a carrier gas also 

limits operation time,^^ although progress has been made by attempts to 

22 ^ 
use H2 or air as a carrier. Recent advances In column composition and 

23 24 

calibration methods have been reviewed by Villalobos and Burchfield. 

Mass spectrometers represent a further advance toward mechanization. 
This method can provide great mass resolution, high sensitivities, or 
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rapid responses, but not all these characteristics can be found in a 

single spectrometer.^^ The necessity of carefully monitoring mass peak 

distributions requires a stable calibration. Analysis may be performed 

real-time, although sensitivity is much greater if preconcentration is 

used. Analysis is becoming more routine as the librax^ of catalogued 

25 

mass spectra grows and better search-and-compare computer programs are 

26 

created to make full use of this collection. Because normal air con- 
stituents are detected along with the trace components, mass fragment 
coincidences with components, such as CO and ethylene, prevent detection 
of these trace species. The efficljency of operation is low because only 
a small part of the information gathered concetms the trace species; this 
must be con^ared to the spectrop ho tome trie techniques, where all the 
collected information is due to the contaminant species alone. Preferential 

diffusion through membranes or ultraflne porous fretted tubes has resulted 

24 

in some sample enrichment. The monitoring of less volatile pollutants 

has been achieved through the use of integrated lon-^ current mass spec- 
27 

trometers. Here again, sample analysis would have to be performed at 

a central location for urban monitoring networks. Extension to operation 

in isolated’ closed environments is hindered by calibration drift problems 

for the more complex mass spectrometers resulting in poor stability. One 

28 

mass spectrometer system for naval applications had to restrict its 
operation to 10 ppm sensitivity for only four or five contaminants to 
meet the stability and dependability requirements. 

The use of a gas chronuitograph as the first stage with mass spectro- 
scopic detection, simplifies the total mass peak intensity pattern and 

15 25 

permits sensitive detection, * but real-time capabilities are lost 
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and the overall complexity and support requirements (weight, power) of 

the system are considerably Increased. One experimental model was 

capable of detecting halocarbons at the 5 parts per trillion level 

29 

with a general accuracy of + 5Z. 

The unlqiie spectral patterns absorbed or emlt^d by each species 

have been utilized as a further step toward completely mechanized 

detection. The only requirement Is the need for frequency discrimination 

In the signal. Over 90Z of all molecules have some absorption between 

30 

2y and 15y In the Infrared. Atoms and homonuclear dlatomlcs, the major 

components of air, do not absorb In this region. Electronic absorption 

bands In the visible or ultraviolet may also be used to identify molecules, 

31 

although the spectra are more complicated and overlaps frequently occur. 

Techniques using grating monochomators for frequency dispersion of 
broadband source radiation may reqiilre non-routine analysis of the absorp- 
tion patterns. In addition, throughput Is small unless resolution Is 
poor. Analysis Is real-time, only seconds being required for a spectral 
scan. If no Interferences between absorbers exist, minimum detectable 

concentrations (MDC's) in the 1-1000 ppm range are possible, depending 

32,33 

on absorption Intensity. * The use of masks at the output of the 

monochromator to check for spectral similarity between the monitored 

signal and the pattern for a species molecule permits analysis, but only 

one or two species may be simultaneously detected. Interference with other 

34 

absorbers is practically eliminated. These correlation spectrometers 
have found wide use, especially In the ultraviolet region, when only a 
few species are to be monitored. These, Instniments, reqiilrlng frequency 
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dispersing elements for use with remote sources or blackbody radiators, 

are subject to mechanical instabilities and are generally too massive 

for spacecraft use. In addition, because of low throughput, sensitivity 
35 

is poor. Good resolution and mechanical stability has been achieved 

36 

by a Fourier transform spectrometer. Single pass £tbsorption yields 

37 

10 ppm sensitivity. These spectrometers combine the properties of 

good throughput, wide spectral range, high resolution, and rapid scan 

times. In addition, signal to noise may be increased by averaging a nundjer 

of scans. Limitations are imposed on sensitivity by the maximum size of 

38 

the mirror and on resolution by the distance the mirror can move. Both 

39,40 36 

urban point sanq>ling * and air-bome Instruments have performed well. 

Use of a multipass White cell permits effective path lengths of several 

39 

hundred meters and a proportional increase in sensitivity. 

Various schemes have been attempted to try to avoid the necessity 
of a monchromator for frequency discrimination. Filters generally have 
too broad a band pass for good discrimination.^^ The iise of a Fabry- 
Perot with high finesse^^ has resulted in very good resolution and it 
is tunable over m of the infrared. NO, NO 2 , CH^, and H 2 O all 

have significant absorption in this region. Gas filter non- dispersive 
infrared techniques (NDIR) have been recently reviewed. In a 
comparison of NDIR with two wet chemical techniques for NO^ detection, 

NDIR had similar sensitivity but was simpler and faster, required much 
less complex analytical techniques, and could give measurements every 2^4 

AC 

seconds.'*'' Although NDIR is widely applicable, mechanized detection is 

limited to 3 or 4 species because a reference cell of each gas to be 

46 

monitored is required. Sensitivities are in the 0.1-100 ppm range. 
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Broadband sources may be eliminated If emission from a reference hot 

gas is used as a species-specific line source. Detection may be by 

monitoring absorption, fluorescence, or thermal heating of the 
44 

gas. The reliability and sensitivity of the best of these systems 

make them good candidates for monitoring the few species unable to 

28b 

be detected by more versatile methods. 

47 

Chemiluminescent monitors also obviate light sources. Their 

^ 48 

extreme sensitivity (ppb. range), rapid response (10 seconds or less) 

are largely counterbalanced by their need for reagents and frequent 

calibration for quantitative detection. Their non-specificity, responding 

49 

to an entire class of reagents, such as NO^, PAN, nitrates, and nitrites, 

also hinders quantitative analysis. In addition, only a single compound 

or class of compounds may be detected in each reaction chainber. Although 

designs for multiple chamber chemiluminescent monitors exist, complexity 

is Increased greatly. It Is unlikely that chemiluminescent methods will 

find wide application. Other technlqties Involving photodissociation 

51 52 

with conductivity detection or magnetic tuning are generally good for 
only a single species and not widely applicable. Lasers suggest them- 
selves as radiation sources because of their extremely high spectral 
53 

brightness. By embodying the frequency discrimination In the source, 
only slnq>le broadband detection is needed. Absorption lines are collisioa 
broadened to about 0.1 cm~^ in an atmosphere of air. Pressure broadening 
and absorption linesha'pes will be considered later 'in the thesis. The 
spectrophotometrlc techniques dlscxissed above were able to achieve resolu- 
tions of 0. 2-0.5 cm ^ at best, and much potentially useful spectroscopic 
Information Is lost. The spectral bandwldths of lasers can be orders of 
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magnltiide narrower than this, and thus they are capable of resolving 

absorption lines of different species, which would otherwise overlap 

and lead to a loss of sensitivity due to interferences. The resolution 

may be further Increased throxigh reduction of sample pressure and 

hence broadening. Gas filter NDIR methods are the only other spectrophoto- 

metrlc technique which similarly Increases in resolution at lower pres- 
40 

sures. Absorption peak intensity is not decreased, and the sharper 
absorption lines result in excellent discrimination. Linewidths con- 
tinue to decrease until the Doppler broadened limit is reached. 

-3-1 

Typical Doppler linewidths are on the order of 10 cm . It should 

be noted that if the molecular line density is such that the lines are 
-3 -1 

less than 10 cm apart, individual lines are not resolvable at any 

44 

pressure and no gain in sensitivity results from pressure reduction. 

Laser Raman methods rely on the inelastic scattering from atmos- 
pheric components being shifted from the laser frequency by an amount 
characteristic of each molecule and proportional to the energy separation 

between the ground and excited states. All atmospheric constituents are 

54 " 

thus detectable by this method. Efficiency of operation is low due to 
the very small scattering cross sections resulting in most of the laser 
light being wasted. Also, a monochromator is needed for discrimination 
of the multi-frequesncy scattered signal. Weight, stability, and effi- 
ciency considerations preclude use of this technique in extraterrestrial 
applications. Sensitivities in the 10-100 ppm range have been achieved 
in point sampling applications.^^ The application of Raman LIDAR for 
remote sensing of pollutants is discussed in a later section of this 
chapter. 
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The obseximtlon of fliiorescence by the gas following an absorption 

of a photon at a particular resonant frequency has also been attempted. 

The reemlsslon of light Is not Instantaneous, but occurs after a mean 

delay, T , called the radiative lifetime of the level. The utilization 

of the large absorption cross sections In the vlslble-ultravlolet and In 

the Infrared make this approach seem promising. Unfortunately, rapid 

Internal energy redistribution occurs among molecular energy levels, 

and the radiated emission occurs over a wide frequency bandwidth, typically 

tens of nanometers In the visible?^ In addition, collisions occurring 

during the lifetime of the excited state may result In loss of excitation 

non-radlatlvely . This process Is called qtienchlng of the excited 

nKilecule. The rate of quenching collisions Is given by 

Q » N Z ) V (1) 

3 2 

where N Is the density of the quenching gas In molecules/cm , O the 

kinetic collision cross section, and v the average molecular velocity. 

Z Is the average number of gas collisions required to quench the excited 

19 

species. In air at standard temperature and pressure (N*2.58 x 10 cm ; , 

9 -1 

the rate of gas-klnetlc collisions (Z = 1) Is approximately 5 x 10 sec . 
The quanttnn yield of fluorescence Is given by 


(^ - ( 1 + Q (2) 

where Is the radiative lifetime. If the quenching rate Is much 

greater than the fluorescence decay rate, very little Isotropic emission 

will be observed. For electronic transitions In the visible part of 

*5 “8 56 58 

the spectrum, x ^ is typically In the 10 to 10 second range. ’ 

Quenching efficiencies have been measured for a number of common pollutanii 

59 

species. Including NO, NH^, and SO 2 . These quantities can show large 
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variations, but most systems of interest have Z's on the order of 1 - 10. 

-5 

For the particular case of optically excited NO 2 , = 4 x 10 sec 

and Z for quenching is approximately 8, so that (j) = 4 x 10 Values 

of quantum yields and effective emission cross-sections for several 

atmospheric contaminants are listed in Table IV. These values are 

for quenching by an atmosphere of ambient air. These cross-sections 
-22 -27 2 -1 

are in the 10 - 10 cm sr range. Fluorescent pollution detection 

is thus several orders of magnitude more efficient than Raman scattering, 
whi^ch has a cross-section of 10 cm^ sr”^ ti^lcally. One elaborate 
scheme using many optical filters was able to detect 0.6 ppb NO^ by its 
fluorescence in the red.^^ 




Table IV 



Fluorescence Quantum Yields for Some Pollutant Molecules 

Molecule 

Excitation Wavelength 

Quantum Yield 

Effective Cross-Section 


in nm 


(Jeff ( cm^ sr“^) 

SO 

300 

1.6 X 10 "^ 

2 X lO"^^ 

NO- 

400-500 

1 X 10"^ 

2 X lO"^^ 

03 ' 

260 

1 X 10 ’^ 

1 X 10 "^^ 

NO 

227 

3 X 10 “^ 

3 X 10 "^^ 

OH 

308 

1.6 X 10 "^ 

2 X 10 “^^ 


For vibrational fluorescence in the infrared, is generally 0.01 

to 1.0 second. Quenching by collislonal deactivation has been measured 

62 63 

for a rnsober of small polyatomic species, ’ for which Z varies between 
2 4 

10 and 10 . Water vapor is notoriously efficient in quenching vibrational 
64 

excitation. A reasonable estimate for Q for Infrared fluorescence 
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would be 5 X 10^ sec and a quantum yield not much greater than 

-5 -5 -4 

2 X 10 should be expected. Values of 10 and 4 x 10 have been 

estimated for NO fluorescence around Its fundamental at 5 
Thus the signal for fluorescent processes may be quite small and laser 
light sources will be required. Power requirements are not as severe as 
for Raman scattering detection, and filters may be tised In place of a 
monochromator. Unfortunately , regions must be f o\md where only, a single 
pollutant absorbs or fluoresces strongly. The necessary diversity of 
laser wavelengths may not be available with a single txmable laser, and 
the desirability of the technique Is decreased, especially for spacecraft 
and submarine operation. In addition, quenching strongly depends on 
the concentrations of various quenchers present and quantitative con- 
centration determination may he difficult. Sample dilution with a gas 
having a small quenching cross-section provides a reduction of the effect 
of other quenchers present. 

The energy not reemltted by fluorescence is degraded to thermal energy 
of the molecules. As we have seen, this may be nearly 100% of the absorbed 
energy. Detection methods measuring the Increased pressure (due to tempera- 
ture rise) In a sample were proposed originally using broadband sources In 
non-dispersive configurations. Use of these spectrophones with l 2 iser light 
sources has resulted In tremendously Increased sensitivity and specificity. 
Detection of concentrations at 0.1-10 ppb level has been achieved using 
high power gas lasers. Early predictions claimed that microphone 

Johnston noise or Brownian motion In the gas cell would be the limiting 

69 

factor and sub-parts per trillion sensitivities would be possible, but 
the presence of atmospheric absorbers, such as water, tends to limit the 
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ultimate sensitivity. Even If no spectral overlap exists, window and 

wall absorption would give rise to a larger signal than Brownian 

72 

motion In the microphone. The ultimate sensitivity of the method 

has been predicted to be 0.1 ppb for air pollutants with strong absorp- 

73 

tlon In a region where water absorption Is weak. Use of an acoustically 

resonant cavity with a feedback loop slaving chopping frequency to a 

72 74 

resonant frequency of the absorption cell has been attempted. * 

12 

Dewey envisions a system which would achieve the same sensitivity as 
Kreuzer,^^ but with the use of much weaker laser sources, thus opening 
the way for the Implementation of tunable sources. The resonant enhance- 
ment results In an Increase In system complexity and a large decrease 
In reliability. Non-resonant schemes must rely on chance overlaps of 
strong pollutant absorption with fixed-frequency hi^ power sources. 

Heterodyne detection has also been proposed^^ as a method for 
Improving detection system sensitivity. A laser is reqiilred to act as 
a stable frequency local oscillator. The signal, thermal emission lines 
from hot pollutant gases, and the local oscillator beam are both 
Incident on an Infrared detector such as copper-doped germanium. The 
detector acts as a mixer and In Its output there exists a frequency com- 
ponent at the difference frequency between the laser and molecular line 
centers. The amplitude of this component Is linearly proportional to the 
signal from the gas.^^ If this intermediate frequency (IP) Is limited 
to frequencies less than 1 GHz by the IF amplifier, a single molecular 
line (at most) may be observed. This IF bandwidth Is on the order of 
the molecxilar llnewldth at atmospheric. If the laser local oscillator 
Is alternately operated at two different frequencies, on and off resonance 
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with the molecxilar species, background oay be eliminated because the 
difference in the signals at the two frequencies will be proportional 
to the molecular concentration. Interference with water vapor emission 
lines and other species would tend to be the limiting factor in detec- 
tion sensitivity. Detection of fluorescence using heterodyne detection 
sxiffers greatly from the narrow IF bandwidth. Only a single line in 

the entire fluorescence band may be observed resulting in utilization 

63 

of only 2% of the fluorescence signal. The noise equivalent power 

of this infrared heterodyne radiometer is given by the formula^^ 

N.E.P. = 2 (hv/n) (Bif/t) ^ (3) 

where v is the detection frequency, n is the quant\aa efficiency of the 

mixer (which may be as high as 0.5), is the IF bandwidth, and t 

is the postdetection integration time. At 10 pm for ■ 1 MHz, 

—16 

T =■ 0.1 seconds, and ri® 0.5, N.E.P. • 2.5 x 10 W. The N.E.P. of the 

heterodyne detection method is actiially greater than for the ideal direct 

detection for equal integration times because of the large value of B^^, 

but in the infrared the N.E.P. for the direct detectors is very much 

greater than ideal due to thermal noise in the detector or circuitry. 

The N.E.P. for heterodyne detection is on the same order as for photon 

counting systems in the ultraviolet. As a result, in the infrared a 

heterodyne detection system may be four to five orders of magnitude more 

78 

sensitive than an optimized direct detection system. Experimental 
investigations^^ were able to achieve detection within an order of 
magnitude of the theoretical limit. Fifty ppm of SO 2 at 390° K was 
detectable in a meter long path. This method has been adapted as 
a remote source monitor for smokestacks and other thermal sources when 
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In the passive mode. Although ambient temperature detection Is possible, 
the method Is not as sensitive as NDIR techniques. Applications using 
a second laser and active monitoring of remote ambient concentrations 
will be discussed In the next section. 

Both opto-acoustlc and heterodyne detection schemes are more sensl> 
tlve than direct detection, but heterodyne methods do not lend themselves 
well to point sampling and opto-acovistlc detection has neither the 

79 

reliability, stability, nor rapid response of the direct detection schemes 
and as a result, the detection scheme chosen would depend on the applica- 
tion. For long duration missions, where reliability Is essential, a direct 
detection system Is obviously preferred. 

Many other proposed detection schemes have severe limitations. 

Microwave rotational spectroscopy requires low pressures to remove pressure 

80 

broadening and large amounts of electrical power. Detection limits are 

In the ppm range. Use of an eibsorptlon cell In the laser results In 

sensitive (sub-ppm) detection, but this process Is severely non-llnear 

8X 

’In absorber concentration and analysis Is difficult. 

82 83 

Magnetic tuning of molecules or of the laser and Stark tuning 
84 

of vinyl chloride has resulted In ppm sensitivities, but these schemes. 

In addition to requiring electromagnets or parallel electric plates , 
are not widely applicable and are generally too complex to warrant their 
Implementation as a single gas detector when NDIR techniques have conquerable 
sensitivity. 

Laser methods offer significant Improvements over blackbody spectro- 
photometrlc techniques for both absorption (because of the additional spec- 
tral Information available) and fluorescence. A direct detection method 
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to measure the differential absorption of laser emission at many frequen- 
cies appears to be very promising both from the standpoint of sensitivity 
and the possibility of application for monitoring of many pollution 
species, t^e will return to considerations of laser detection capabilities 
when we dlsctiss remote monitoring applications below. 

The final point sampling technique that will be considered is a hybrid 
method involving spectrophotometry and mass spectrometry. Most mass spec- 
trometers rely on an electron beam colliding with the sample and ionizing 
it. The fragments are detected by their mass/charge ratio after suitable 
focusing by magnetic or electric fields. If, however, photoionization 
is used, frequencies may be selected which do not result in the ionization 
of the major atmospheric constituents. Such a scheme has obtained ppm 

sensitivities, and envisions ppb MDC levels with ionization frequency 
85 

maximization. 

The use of two-step photoionization involving a tunable Infrared 

photon and an ultraviolet photon not individually capable of ionization 

86 

has been suggested. Tuning of the infrared frequency would result in 
changing mass peak patterns in the mass spectrometer output. This two- 
dimensional variability could be used to distinguish substances with 
nearly identical mass peak patterns, but different Infrared absorption 
spectra. 

D. Remote Monitoring Techniques 

Remote monitoring schemes must be based on spectrophotometric 
methods. Extension of these methods, however, is not straightforward. 
Broadband Incoherent sources soon meet Intensity limitations. Substitution 
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of solar or sky backgromd radiation as a light source has had mixed 

V 

87—91 

success. Unfortunately, with these light sources, operation is 

restricted to daylight hours, and performance is severely hampered by 

87 

rapid changes in sky brightness. In addition, dispersive instrum en ts, 

with their low throughput, are confined to low resolution work to gain 

88 

sufficient signal intensity. If direct solar radiation is used, the 

89 

system's spatial variability is destroyed. Nonetheless, sensitivities 

89—91 

in the parts per billion range are achieved for many substances, 

and light source Improvement is the obvious remedial action. Here, as 

in the point sampling applications, the narrow frequency bandwidth, high 

intensity, and easy beam colllmation recommend the use of laser light 

sources. All the proposed monitoring systems rely on atmospheric 

attenuation of the beam giving rise to frequency discrimination. Atmos- 

92 

pheric attenuation may be singly described by the Lambert relation 

I/^o " ®*P {- /Qa(x) dx } (4) 

where a(x) is the local attenuation coefficient per unit length and 

I/I is the fraction of radiation remaining in the beam at distance L. 
o 

Atmospheric attenuation has three major sources: Elastic scattering by 

atmospheric gases, elastic scattering by aerosols, and absorption by 

gaseous constituents in beam path. Inelastic Raman scattering by 

molecules gives rise to a small attenuation also. 

The atmospheric elastic scattering component due to scatterers 

much smaller than the wavelength of the incident radiation is called 

Rayleigh scattering. It is primarily attributed to scattering by 

93 

atmospheric gaseous constituents, and is given by the formula: 
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o o e 

where u Is the center frequency of the nearest resonance, q Is the 
0 

charge of the electron, Is Its mass, and Is the free space permit- 
tivity. Rayleigh scattering is most lsq)ortant in the ultraviolet due to 
Its Inverse A dependence. Values for the Rayleigh elastic scattering 
coefficients, different wavelengths are given in Table V. 



Table 7 

Rayleigh and Mle Scattering 

Coefficients 


X 

(micron) 

Og(X) 

(km 

Oj|(X) 
V-1 km 

(km“l) 

V“5 km 

V=>10 km 

0.2 

0.70 

7.0 

2.2 

1.4 

0.3 

0.14 

5.6 

1.4 

0.84 

0.4 

0.044 

4.7 

1.1 

0.58 

0.6 

0.0083 

3.7 

0.72 

0.35 

0.8 

0.0026 

3.1 

0.54 

0.24 

1.0 

0.0006 

2.8 

0.43 

0.18 

5.0 

0.000001 

1.1 

0.086 

0.024 


0.00000006 

0.72 

0.043 

0.010 


The elastic scattering by particles or droplets of size comparable 


to the wavelength of light Is called Mle scattering. Because a cotq>lex 


size distribution exists for particles between 0.1 and 10 Um, use is 


made of an approximate formula relating the Mle elastic scattering 

94 


coefficient, oi jj (X), to the visual range' 


a ^(X) 


3.91 


0.55 


0.585 V 


1/3 


km 


-1 


( 6 ) 
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V is the vlsijal range in kilooeters and X is the wavelength In microns. 
Hie scattering is much less wavelength dependent than Rayleigh scattering. 
In contrast, Mle scattering is much more a function of changing atmos- 
pheric-conditions than Rayleigh scattering. Mie elastic scattering 

coefficients, a (X), are also given in Table V for different X and V 
M 

values. The combined atmospheric elastic scattering coefficient, a g(X), 

ct^( X) « X) + X) (7) 

is plotted as a function of wavelength for various visibilities in 
95 

Figure 1. 

Atmospheric attenuation due to absorption of radiation is related to 

92 

the concentration of the absorber present by Beer's Law 

— - exp /q<J(^) N dx » exp ^-/qE(X)c dx (8) 

where is the fraction of light remaining at distance L, CT ( X ) is 

2 

the absorption cross-section (cm ) , N is the number of absorbers present 
per cm , £ ( ^ ) is the molar absorption coefficient per unit concen- 

tration per unit length, and c is the molar concentration. The elastic 
scattering processes have a much more slowly varying frequency dependence 
than the absorption cross-sections, and thus it is possible to txme 
across several absorption features maintaining the same elastic back- 
gro\ind attenuation. However, overlap may exist emd several species may 
absorb radiation of the same wavelength. Infrared monitoring is limited 
to spectral regions where absorption due to the natural atmospheric 
constituents CO^ and H^O is small. Below 250 nm in the ultraviolet O2 

absorption becomes severe and limits operation distances. Atmospheric 

44 58 

constituent absorption Interferences are shown in Figure 2. * 
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Figure 1; Atanospheric Scattering as a Function 
of Wavelength and Meteorological Visibility for a 
Number of Different Wavelengths (from Ref. 95) 

Figure 2j Atmospheric Constituent Interferences 
and Atmospheric Transmission Windows (from Ref. 58) 

Figure 3^; Background Sky Radiance as a Function 
of Wavelength (taken from Reference 58) 
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FIGURE 





FIGURE 2 



PAN 





Both single and double-ended absorption monitors have been proposed. 
Single^ended operation relying on remote reflectors is preferred because 
only a single detection system is required. Spatial inhomogeneities are 
not detected, only the average concentration in the absorption path is 
monitored. Thus we have for this monitoring system 


96-110 


I . *• I , £ exp 

sig eff ^ 


-2LY a.(X) N. 

• i i i 


^eff 


|-2La^(X) ' 


with 


eff 


I f exp -2 a ( X )L 


(9) 


( 10 ) 


where is the measured intensity of the returned light, L is the 

— -3 

distance to the reflector, is the average concentration (cm ) of 

absorber 1, a ^( X) is the wavelength dependent absorption cross-section 

for species i, accounts for the reduction of signal due to atmospheric 

elastic scattering and loss at the remote reflector, f being the fraction 

I 

of incident intensity returned by the reflector, and C is a collection 

efficiency factor dependent on beam spread and detector size and efficiency. 

The wavelength dependence of cr( X) will be different for each absorber, 

and by tuning across a sufficLeEtly wide wavelength region, individual species 

would be identifiable. Proper choice of wavelength would minimize or remove 

absorption overlap simplifying detection of that species. Placement or 

choice of reflectors at radially increasing distances from the monitoring 

96 

station would permit spatial resolution. Objects such as trees and 
97 

buildings have been used as remote reflectors for laser beams. Reflec- 
tivities of over ten percent have been achieved using these non-cooperative' 

111 -17 -21 

reflectors. Because absorption cross-sections are large, 10 - 10 

2 112 

cm , this method should have good sensitivity even for low laser powers. 
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Before an evaluation of the limitations of this method, laser backscatter 
techniques will be considered. 

Monitoring schemes which rely on atmospheric attenuation to provide 
Intensity which Is scattered back In the direction of the detector have 
been widely proposed. For this situation, the measured Intensity of the 


returned light, is given by 


slg 


leff 


(- ♦ “ a <»1 •• " 


where now 




-{Oj(X) + aJX)} l1 


out 


( 11 ) 


( 12 ) 


I J.J. Is the total Intensity of the backscattered radiation, and Is 
eff 

dependent on the amount of Intensity reaching the sampled region 


do exp 


( )’> 


on the length of the sanq)led region (D) and on the 


magnitude of the scattering process ( a ^( X ) ). Once the radiation 

Is scattered, further atmospheric attenuation, both elastic and absorp- 
tion, occurs between the scattering region and the detector. Because 
the radiation is scattered Isotropically, the factor C is introduced, 
which takes Into account the solid angle subtended by the detector and the 
detector efficiency. Range resolution is possible by gating the detection 
electronics to monitor the backscattered signal at varying delays 
after the laser pulse. Three different LIDAR approaches to res»3te 
monitoring exist; Raman, fluorescence, and differential absorption. 

The Raman Inelastic scattering cross-section is _exl^remely small, 

3 

typically 10 less than Rayleigh cross-sections, and thus Raman processes 
were not considered to give rise to significant attentuatlon of the 
beam. Like Rayleigh scattering, Raman backs nattering has an inverse 
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where v Is the wavenumber of the exciting beam, Is the ILaman wavenumber 

shift, g Is the degeneracy, h and k are the Planck and Boltzmann constants, 

t ' ’ 

T Is the temperature of the molecules, a is the Isotropic part of the 
derivative of the polarizability tensor and y is the anisotropic part of 
that derivative. P is the Incident laser power.. 

Thus , 

Ot.(X) - cr^(X) N (14) 

3 * 

— —3 

in equation (12), where N is the average concentration (cm ) In the sampled 
region. The X value used In equation (12) is the laser wavelength. The 
wavelength . of .the backscattered iradiation ?ln .equation .('11>) .is tthe ishlited 
wavelength, the shift being characteristic of each atmospheric constituent 
present. Some molecular absorption, a ^ > of backscattered 

radiation may occur in addition to atmospheric elastic scattering, further 
complicating determination of pollutant concentrations . Use of a mono- 
chromator or filters at the central detection location would allow possible 

54 

detection of all trace species. Because tunabillty Is not necessary, 
more powerful fixed frequency lasers can be used. The v** dependence 
favors work in the ultraviolet. If an electronic transition of one of 
the gaseous pollutants happens to be nearly resonant with the laser 
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frequency, an enhancement of the scattering cross-section of several 

119 120 

orders of magnitude can occur. ’ The Raman cross-sections of 

several gases at different wavelengths are given In Table VI. Resonant 

Raman cross-sections are also shown to demonstrate the enhancement. 

Unfortunately, the resonant enhancement would be In effect for only, at 

most, a few species for any given wavelength, and a ttmable source, 

which has lower power levels, would be needed. Because Raman processes 
121 

are instantaneous, the laser pulse width alone determines the limit 

of depth resolution by this method. The overlapping of the Raman 

spectra from various molecules results in loss of some sensitivity. Use 

of the rotational Raman effect, which has a cross-section a factor of 100 

larger than the vibrational Raman effect, has been suggested as possible 

through use of a Fabry-Perot filter to eliminate pump laser light and much 

of the background. No monochromator would be needed. The free spectral 

range of the etalon would be the rotational line separation. A different 

122 

etalon would be required for each species to be detected. 

The remote monitoring of resonant fluorescence emission has also 
been suggested. Here H^O and particulate quenchers cannot be removed, 
and in addition, the fluorescent emission bandwidth requires a wide 
receiver bandwidth, making background intensity more of a consideration. 
For fluorescence monitoring (A) equation (12) is 

( A ) » a (X ) (j> N (j^5) 

where Cj(X) is the species absorption cross-section at the laser 
wavelength, (p is the quantum yield defined in eqtiation (2) and N is the 

3 

average density of the monitored species in molecules/cm . For this 
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Table VI 


Raman and Resonance Raman Cross-sections 


Molecule 

Raman Cross-section Wavelength Ref. 

Res. Raman Cross- 

U9„ , 

Wavelength 


( 2 -Iv 

(cm sr ) 

(nm) 


section (cm^ sr 

(na) 


2.5 X 10“^° 

347.1 

123 



z 

3.0 X 10"^° 

347.1 

124 




2 X 10”^^ 

694.3 

115 



N 

3.0 X 10"^° 

347.1 

123 



z 

2.5 X lO"^ 

347.1 

124 




1.9 X 10"^^ 

694.3 

115 



HO 

9.5 X 10“^° 

347.1 

123 



z 

6-8 X 10“^° 

347.1 

124 



N02 




5.6 X 10"^^ 

454.7 

NO 

-30 

0. 2-1.0 X 10 ^ 

347.1 

124 

-23 

1.3 X 10 

226.2 


-29 

1-5.5 X 10 

337.1 

56 



CO 




/ c in“30 

4.5 X 10 

206.3 

SO, 

-29 

1.3 X 10 ^ 

347.1 

124 

1.0 X 10‘^^ 

300.0 

z 

-31 

3.0 X 10 

694.3 

115 




-29 

4.8 X 10 ^ 

337.1 

56 



NH^ 

-29 

1.0 X 10 ^ 

347.1 

124 




-29 

1.4 X 10 

347.1 

124 



CO, 

-29 

5.5 X 10 

337.1 

56 



Z 

” -30 

3-4.5 X 10 

347.1 

124 




7.3 X 10“^^ 

694.3 

115 



N2O 

7.3 X 10"^® 

347.1 

124 
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method A In equation (12) is the laser wavelength and X In equation (11) 

Is the wavelength band of the output fluorescence. Molecular absorption 
In equation (11) by species present. Including the fluorescing one, could 
give rise to too low remote concentration estimates. Elastic scattering 
losses also reduce the returned signal. Laser wavelength constraints are 
the same as for absorption, a frequency must be fotmd which has a signifi- 
cant absorption cross-section. Excited state lifetimes, even with the 

64 

presence of quenching collisions, are several microseconds In the Infrared, 
limiting range resolution to 1 km or more. In the ultraviolet, lifetimes 
are shorter and resolution may be limited to laser pulse duration or 
electronics gating Interval. 

The differential absorption technique relies on atmospheric elastic 
scattering to provide a remote light source which Is relatively constant 
over a range of wavelengths. The differential molecular absorption of this 
retximed signal Is then used to determine the average concentration of 
the absorbing species between the observed scattered region and the detector. 
Thus 

a, (X) = a_ (X) (16) 

j =• 

and the laser and backscattered Intensities are at the same wavelength. 
Elastic scattering occurs Instantaneously and the depth resolution Is limited 
only by the laser pulse width. By varying the gating delay, different 
atmospheric regions are used as the remote source and a radial map of 
absorber concentrations can be constructed. Here again a tunable source 
Is required and monitoring must be performed In a region where there Is 
a structured absorption pattern. The sensitivity of the method Is Increased 
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as the difference in absorption crass-sections between two nearby wave- 
lengths increases. Substituting in equations (11) and (12) for two 
wavelengths X and X* and ratioing the results, we see 




(- 


eaq) -2 (o - a ) N L 


(17) 


if only a single species absorbs ^ and X ' . Even if several species 
absorb, monitoring of a sufficient number of frequencies should allow 
determination of all absorbers present. The value of ®ay be 

small on clear days or at larger distances on foggy days, resulting in 
a low signal level returned to the detector. As in the absorption , 
method, the returned signal is at the laser frequency and narrow band- 
width detection is possible. 

We now have a sufficient basis for comparing the remote detection 
schemes. M^y conq>arlsons have been performed previously, 
but none considered all of these methods together or attempted to conq>are 
them with as few variables as here. Rewriting the equations for each 
of the proposed monitoring schemes yields 


slg 


Iq exp{-2 a^L} f C exp {-2a^L} 


E 

^e' 


Absorption 'With reflector 


’■sig “ ^o |f C'(.05) exp{-2 a^L} 


Ig^g “ Iq exp|-2agL} D C exp{-a^ l} a(X) (J)N 


xp{-2 ttgL} DC aj^(X) iij 


out 


Absorption with non- 

cooperative target 

Fluorescence 
backscatter 


I , « I e 
sig o 

I . » I exn(-2a,Ll fp 
sxg o - " ^ ' i. 


Henceforth, we let S 


C exo(-2 a lV a. (X') 
A ' E 


Raman backscatter 

Differential absorption 
backscatter 


4 
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The asstimptlons that have been made are: 




“a ’out * ’back 

“A'^’back-O 


for Raman and fluorescence 
backscatter. 

0 for Raman backscatter. 

for fluorescence backscatter. 


4. A 5% reflectivity from the non-cooperatiro, .remote reflector. 

Assumptions 2 and 3 would give rise to a larger return signal than would 

be the case If atmospheric absorption was not zero. The factor 

exp {”2 a £ (X)L} is common to all schemes and will be neglected 

teiiq>orarlly. Let us now determine the relative returned signal from 

a pollutant present at low concentrations, 1 km distant. Range resolution 

will only be 300m, finer resolution will result in proportionally 

“■28 2 

decreased signals. The value of a (X ) Is chosen as 10 cm ?■ 

R 

""I —19 2 58 

molecule steradian ^ o (X) as 3 x 10 cm / molecule steradlan, 

—2 —6 

and as 10 to 10 . Values of the terms compared, S, are given In 

Table VII for 1 ppm and 10 ppb average pollution concentrations, and for 

different remote distances from the detector. Use of 10cm radius 

collection optics and 10% photon to signal efficiency gives 
-9 ' 

C s 1 X 10 at 1 km. fC has been previously estimated as less than 
-4 

10 , depending on beam colllmatlon and reflector size, using a 


similar scheme to that above.' 


As expected, the return signal for absorption Is many orders of 
magnitude greater than the rest. For all the absorption schemes, not 
the total signal returned, but the change In signal as a function 
of wavelength Is lnq)ortant. Thus, If a 1% change In signal Is Just 
barely detectable, and the change in absorption coefficient in 
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Table VII 


Comparison of Remote Monitoring Techniques 

Method Concentration Values of S for Remote Ranges, L, of . 

- 10 km 1 km 100 m 10 ai 


Absorption 

1 ppm 

-12 

1.7 xlO ^ 

2.1 xl0“^ 

8.6 xlO"'^ 

9.8 xlO 

with Reflector 

10 ppb 

8.6 xlO**® 

9.8 xl0“^ 

1.0 xlO"^ 

1.0 xlO' 


Absorption with 

Topographical 

Reflector 


1 ppm 
10 ppb 


8.5 

xlO"^^ 

1.1 

xlO ^ 

4.3 xlO”^ 

4.9 

xlO' 

4.3 

xl0~^ 

4.9 

xlO"^ 

5.0 xlO"^ 

5.0 

xlO 


Raman 1 ppm 7.8 xl0“^^ 7.8 xl0"^° 2.6 xl0“^® 2.6 x10 

Backscatter 10 ppb 7.8 xlO"”^^ 7.8 xlO 2.6 xlO 2.6 xlO 


Fluorescent 

Backscatter 


1 ppm <j)«'10 

-2 

1 ppm (fja'lO 
lOppb ^**10 ^ 
lOppb (^s»10 ^ 


9.7 


1.1 xlO”^® 

2.2 

X10-1‘ 

-12 

2.3 xlO 

9.7 

.10-1* 

-12 

1.1 xlO 

2.2 

zao"^“ 

2.3 xlO"^ 

2.2 

:ao-2« 

2.3 xl0“^® 

2.4 

*io-« 

2.4 xl0“^^ 

2.2 

id.O''-* 

2.3 xlO"^^ 

2.4 

tL 

o 

1 

I-* 

to 

2.4 xl0“^° 


Differential 

1 ppm 

5.0xl0"^\ 

6.3x10T^^V 

8.7xl0”^V 

9 . 7xl0"^' 

Absorption 

10 ppb 

2.6xl0"^V 

2.9xlO"^°V 

1.0xl0"®V 

l.OxlO"^ 


of Backscatter where V =• 5.1 for visibility of 1 km at 400 nm. 

V = 1.5 for visibility of 5 km at 400 nm. 

V - 0.043 for visibility of 5 km at 10 ym. 
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-19 2 

equation (17) Is 1;5 x 10 cm /molecule, steradlan, we have 


I , (X^) 

slg 




exp{-2( a'(X) - a (X)) NL] 


min 


ln(0.99) “ -2(1.5 x 10 ) cm /molec. *str (NL) 

- 3 

(NL) ^ » 1.2 X 10 ppm-cm 


( 17 ) 


(18) 


and If a mlnlmim range resolution of 300m Is desired, 0.O41 ppm Is the 
minimum detectable concentration for this resolution for any distance 
where Ig^g ( X ) la significantly above background. 

As a comparison, the standard wet chemical Saltzman method would require 

six hours of preconcentration In order to achieve this sensitivity. The 

sensitivities of various methods for detection of NO^ are presented In 

Table VIII. The remote methods are all listed for 300m range resolution. 

It Is seen that the differential absorption method provides comparable 

or superior sensitivity to the other remote monitoring schemes, and Is 

capable of detecting ambient NO 2 at a level below the EPA Maximum Exposure 

Level for this resolution. Point sampling methods are currently more 

sensitive, but If remote monitors are capable of sufficient sensitivity 

to monitor ambient concentrations, a single remote probe will obviate an 

entire network of local point samplers. Currently, equivalency standards 

130 

between path and point monitors are being set. Once these relations 
are established, remote monitors may be used to establish the adequacy of 
a point monitor to represent the ambient conditions In Its local area. 

We will now consider the range capabilities of the various remote detection 
schemes. 

Photon limited detectors In the visible have minimum detectable powers 
of 9 X 10 W, If background radiation power Is comparable or smaller In 
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Table VIII 


Detection Capabilities for NO^ for Various Techniqijes 


Method 

Detection Limit in ppm 

Reference 

Point Sampling: 

Saltzman (colorimetric) Method 

0.24 

45 ,127,129 

with 1 hour preconcentration 

Phenoldisulfonic Acid 

10 

128 

Derivative Absorption Spectros- 

0.01 

33 

copy with multipass cell 



3 

Chemiluminescence 

0.01 

48 

Fluorescence 

0.001 

61 

Optoacoustic 

0.0001 

67 

Remote Monitoring: 


- 

Non-dispersive Infrared 

10 

90b 

Fluorescence 

1 

60 

Differential Absorption of Back- 


0.15 

145 

scattered Radiation (visible) 
Differential Absorption of Back- 


0.2 

78 

scattered Radiation with Hetero- 

dyne Detection (Infrared) 



Absorption 

0.04 

Predicted, 
this work. 
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( 19 ) 


53 

magnitude. ^ee Appendix ij The background power received at the 
detector for an extended source filling the receiver field of view is 

Wkground " a^(X))L} AX V A N(X) 

where AX is the optical bandwidth of the system, V Is the field of 

view of the collection optics, and N(X-) is the "Spec-t-ral radiance of 

the background source. The spectral radiance of the sky under clear 

S3 

daytime conditions is shown in Figure 3. In the visible 
_2 —2 —1 —1 

N ( X )° 10 watts cm ]im steradlan . For the absorption schemes 

3 o 

the optical bandpass may be very narrow (10 ym) and for 1° field of 


1.5 X lO"^^ W 


view and for a 1 ysecond gate time (300m resolution), 
background 

for 5cm radius collection optics, and these methods are not background 
limited. Assuming a 100 KW, 1 ysecond tunable laser in a wavelength 
region where significant absorption changes occur we have, substituting 
in equations (11) and (12) using Table VII, for absorption using the 
atmospheric attenuation coefficients in Table V 

Ig^g (X) = 10^ Wx4.6xl0 ^x2.1xl0 ^«9.7 xlO b for 1 ppm V=5km 

“4.5 X 10 b for lOppb V“5km 

for Ikm.I^^^(X) values for absorption, and differential absorption of 

backscatterlng are given in Table IX for various distances and visibilities. 

Note that I . = I, , j at much smaller distances for the backscattered 

sig njackground 

signal. Note also that both methods are useful at greater distances in 

the infrared. This calculation does not take into account atmospheric 

absorption by atmospheric gases in the infrared where total atmospheric 

131-133 

attenuation is about 1% per kilometer and poses no serious 
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Table IX 


Comparison of Absorption and Differential Backscatter Methods 


Method Visibility Wavelength Signal in Watts for lOppb Pollutant 

(km) (Jim) Concentration over Ranges of 

10 km 1 km 


Absorption 

10 

0.4 

• 

00 

X 

Os 

1 

o 

1.3 



with 

5 

0.4 

0 



4.5 

X 

lO”! 

Reflectors 

1 

0.4 

0 



3.4 

X 

lO'"^ 


10 

10 

7.0 

X 

10"^ 

9.6 




5 

10 

3.6 

X 

10"^ 

8.9 




1 

10 

4.8 

X 

10~^ 

2.3 



Differential 

10 

0.4 

5.4 

X 

10-16 

3.9 

X 

io-« 

Absorption of 

5 

0.4 

0 



2.0 

X 

10-‘ 

Backscatter 

1 

0.4 

0 



5.1 

X 

io“^ 


10 

10 

2.1 

X 

10“^ 

2.8 

X 

lo”^ 


5 

10 

4.7 

X 

10”® 

1.1 

X 

10-® 


1 

10 

1.0 

X 

10“1^ 

4.9 

X 

VO 

1 

o 


limitation to beam transmission. 

Using a similar laser to 'detect fluorescence backscatterlng xylelds 
5 -2 -1^ 

I . (X ) *• 10 W X 4. 6 X 10 X 1.1 X 10 * 5.1 X 10 for Ippm V=5km 

=H 

» 5.1 X 10 for lOppb V=»5km 

However, for fluorescent reemlsslon a broader observation bandwidth is neces- 

64 

sary, unless much of the signal Is to be lost. Using a 50nm bandwidth 
in equation (18) yields background ” * lO”^^ W/ Jisecond. Limitation 

of the gate width to 1 jisecond may also result in the loss of some signal 
Intensity. It is seen that although a single pulse may be able to detect 
1 ppm (S/N a 1), 10^ pulses, at least, and signal averaging are needed for 
10 ppb sensitivity, and real-time aspect may be lost. 
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Raman backscatter may make use of a fixed frequency laser having a 
10 MW, 100 ysecond pxilse. Then 

I , (X)= (10^W)x4.6xl0"^x7.8xl0“^° = 3.6xl0“^V for 1 ppm V = 5km 
sxg 

Increasing size of the collection optics will result in an Increase in 

background radiation also. By Increasing collection optics size to 

o 

50 cm radius, while decreasing field of view to 0.01 (2 x 10 


sig 


2 3 

Is Increased by 10 . Nevertheless, 10 pulses are required to 

.7 


steradlan) the background radiation power may be held constant while 
I 

achieve 1 ppm sensitivity, 10' pulses for 10 ppb sensitivity. All real- 
time Information would be lost in the time required for this large 
number of pulses of these low repetition rate lasers. It, has. In 
addition, been assumed that all background Interference by elastic 
atmospheric scattering was removed by the monochromator or filter. 

3 

Because elastic scattering cross-sections are 10 larger than Raman cross- 
sections, and the Raman signal arises from only 10 ^ of the moleciiles 

9 

present (1 ppm), the elastic signal may be 10 larger and not be com- 

58 

pletely removed by frequency filtering. The resulting sensitivity 

Increase for resonance Raman scattering Is at least partially canceled 
2 

by 10 lower assumed power for a tunable laser. Raman schemes have had 

some success at detecting remote sources of pollution up to a kilometer 

distant with sensitivities In the 0.1 - 100 ppm range depending on the 

134-139 

complexity of optical detection system. 

A few fluorescence backscatterlng schemes have been proposed,^^’^^*^^ 
but none as yet has been Implemented. Differential absorption techniques 
have been optimistically touted as a method providing range resolution 
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Several experimental Investi- 


_,*u 1 * ^ 1 1 1 97,98,1A0-143 

with large return signal levels. 

gatlons using tunable visible lasers have resulted in 0.1 - 1 ppm detection 

144 145-147 

limits for 300m resolution for SO^ and 0^ and NO 2 . Alternate 

operation at two wavelengths was used for differential absorption. Simul- 
taneous operation on two different wavelengths 1.3 nm apart has been 

. achieved through use of an intra-cavity Glan-Thompson prism and two different 
97 

gratings. In this configuration, atmospheric turbulence effects, as well 
as attenuation effects, are minimized. Nevertheless, because of the rela- 
tively small backscattered signal, workers using this method required 11 

hours to map the concentration of NO^ In a 300m x 400m area, 750m distant, 

147 4 

with a 100m resolution. Over 4 x 10 laser pulses were required. Obvi- 
ously, mapping time could be decreased by technological Improvements giving 
rise to a faster pulse repetition rate, but this method currently seems far 
removed from a city-wide monitoring network using a visible laser. With 
the application of heterodyne detection techniques, detection quantum effi- 
ciencies in the infrared approach those of visible detectors. Here again 
mlnlmim detectable concentrations of pollutants are expected to be in the 
0.1 ppm range. In the Infrared interference from water vapor lines 

tends to limit the sensitivity of the method since atmospheric scintilla- 

149 

tion is much less than in the visible. In view of thls^ several authors 

have claimed ultimate sensitivities in the 10 ppb range for regions free 

, . ^ X 150,151 

of interference. 

Use of heterodyne detection for passive remote monitoring has been pro- 
152 153 

posed * and results in a sensitivity of 0.1 - 10 ppm. A similar 
Increase in range and sensitivity would result if heterodyne detection were 
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used in a differential absorption monitor. The additional espense and 
Increase In complexity Is probably not wairranted for this technique, which 
already has excellent capabilities In both these areaa. 

E. Conclusions 

It has been demonstrated that absorption methods are capable of sensi- 
tivities comparable or greater than other remote monitoring schemes, while 

having a much larger signal returned to the detector. Sensitivity Is 

39 44 

Increased with distance sampled. The use of multipass cells * should 

permit similar sensitivities for point sampling applications. In this 

situation where background may be neglected, low power cw lasers may be 

used. The low power requirements of the absoirptlon techniques make it 

especially desirable. Eye safety maximum permissible exposure (MPE) set 

-7 “2 

by the American National Standards Institute are 5 x 10 J cm In the 

jg 62 

visible, and 0.56c J cm in the Infrared. The MPE value in the 
visible will severely restrict use of high energy pulse lasers In urban 
applications. The IR standard may be easily met for both pulsed and cw 
operation. 

Because the absorption technique appears promising both for point 
remote monitoring, and because the minimum absorption requirements are ruvt 
stringent, there exists the possibility of being able to use a single 
laser source to simultaneously detect many atmospheric contaminants. Experi- 
mental work has been undertaken in order to further evaluate the deslraaiiity 
of this method. 
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II MONITORING COMPLEX MIXTURES OF TRACE CONTAMINANTS 

A. Introduction 

It was seen in the previous chapter that a laser nionitoring scheme 
which Is capable of detecting a single substance with no Interference' 
from other pollutants present may be desirable as a species-specific 
detector, but such a scheme would require a complex monitoring network 
of many different lasers to adequately nmnitor ambient air quality. 

A necessary consequence of the versatility to be gained by operation in 
the infrared is the possibility of spectral overlap between different 
absorbing species. Due to spectral complexity, a reference frequency 
marker is needed to insure reproducibility and accuracy if a continuously 
tunable laser source is used. All other frequency measurements could 
then be made with reference to this frequency. On the other hand, if a 
line-tunable CO2 laser is used, the laser lines may be carefully stabilized 
to the zero velocity molecular transition frequency, and reproducibility 
ensured. Knowledge of the absorption coefficients of each gas at 
selected laser frequencies will permit the absorption of a mixture of 
these gases to be assigned to the relative abundance of each absorber, 
provided a sufficient number of frequencies are used. 

This project was to determine the feasibility of using a frequency 
stabilized CO2 l£iser for the detection of each trace component in a 
mixture of absorbers and to find the number of frequencies required for 
sensitive reproducible detection. The experimental work consisted of 
three parts: 

1 . Frequency stabilization of the CO2 flowing gas laser; 
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2. Setting up a detection system capable of detecting a wide 
range of attenuations to permit accurate absorption coefficient determina- 
tion for a nundber of likely contaminants ; 

3. Analysis of mixtures of these gases to test the detection 
capabilities of this absorption method. 

A CO 2 laser was chosen becaxise the laser frequencies, which span the 

929 cm ^ to 1082 cm”^ region of the infrared, are in an atmospheric 

transmission window. Water vapor absorption Is Insignificant for short 
131 

path lengths. Molecular oxygen and nitrogen do not absorb anywhere 
in the Infrared. Thus, absorption will be due solely to atmospheric con- 
taminants, providing high Information content In the detected signal. 

217 

Because over 100 laser transitions have been observed for the CO^ laser, 
there is a large degree of freedom in number and choice of wavelengths 
to be used for mixture monitoring. 

B. Laser and Detection System 

The CO 2 laser gain medium was a water-cooled low pressure flowing 
gas electrical discharge sustained by a current regulated 10 kV power 
supply. The laser cavity is composed of a curved mirror (f = 5m) with 
a centered output hole and a grating blazed for lOym (Bausch aind Lomb). 

An intracavity iris diaphragm near the front mirror confines lasing to 
the TEM^^ mode on a single transition and controls output power levels. 
Instabilities in the laser frequency and intensity arise from linear 
expansion of the aluminum support girder due to thermal fluctuations. 

As a result, the frequency of the cavity modes are shifted. The 
resultant laser frequency is pulled off the peak of the gain curve 
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to a freqtiency Intermediate to the frequencies of peak gain In the laser 
medium and cavity gain. To Insure that the CO^ laser frequency is stable, 
the maximum of a cavity mode is kept centered on the gain distribution 
of the medium. This is done by mounting the cavity mirror on a 
piezoelectric transducer. A reference sine wave is generated by a 
lock-in amplifier (PAR Model 120) , is then amplified by a KEPCO OPS 2000 
Op-Amp and applied to the transducer. The resultant mechanical oscilla- 
tions of the mirror produce synchronous fluctuations in the laser 
intensity as the cavity mode, (and hence the laser frequency) is swept 
across a small part of the medium's gain curve. The laser intensity 
fluctuations are monitored by an InSb detector (Opto-electronics Mullard 
type ORP-10) looking at the laser output off the normal to the grating. 

The AC part of the detector signal is proportional to the slope of the 
gain curve and is smallest when the oscillation of the mirror is cen- 
tered on the maximum of the gain curve, which corresponds to the zero 
velocity molecular transition frequency. The magnitude of this AC 
signal is the output of the lock-in aa 5 >lifier when the detector signal 
is the input. The output of the lock-in is Introduced in the high 
voltage KEPCO Op-Amp to produce a bias voltage on the transducer suffi- 
cient to center the mirror on the maximum of the gain curve. 

218 

This circuit was originally Implemented by Nowak, but it was 
discovered that the Spectra-Physics 510 Transducer used in his set-up 
was designed for positive voltage. The KEPCO supply is capable of 
producing negative high voltage only. It is doubtful that the crystal 
ever responded correctly to voltage of the wrong polarity. A new 
transducer (Burleigh Model PZ-90) which was conq>atlble with the supply 
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was ordered. It was then discovered that the supply was capable of 


outputting up to a 3 KV transient when turned off or to 2 K7 If driven 
by sufficient Input. The transducer was designed for 1250 7 maxinnni, 
either polairlty. Thus It was necessary to set up a network which would, 
limit the maxlmxnh output voltage and allow the contribution from each of 
the inputs to the Op-Amp to be of the correct magnitude. The circuit Is 
shown In Figure 11. The voltage divider at the output prohibits the 
voltage seen by the transducer from exceeding 1130 V for a 2000 V maxi- 
mum output. These resistors permitted, at most, 2 ma to flow out of 

Op-Amp and do not seriously load the Kepco amplifier or hamper Its 

frequency response. The gain of the Op-Amp for each of the Inputs 

219 

could be calculated using the relation 

G - (22) 

where R^ is the total resistance at the output of the Op-Amp, and R^^ 

Is the effective input series resistance. The response of the transducer 

was 8 |jm/1000 V. AS ym mirror translation would result in an Integral 

218 

change In the number of wavelengths in the cavity. Allowing the lock-in 
output to cause, at most, a shift of half this amount would- permit maximum 
compensation capabilities while ensuring the two ca'vlty modes woiild not 
be simultaneously competing. A series resistance of 16.56 was intro- 
duced Into this Input to the Op-Amp and correspondingly the maxi mum high 
voltage output that could arise from the lock-in signal was 310 V (2.5 ym 
translation) . The sine-wave reference signal was amplified to give rise 
to a 0.5^m peak to peak oscillation of the mirror at the reference fre- 
quency. The magnitude of the amplified signal from the manual voltage 
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Figure 11 ; Frequency Stabilization Network 

Figure 12 ; Lock-In Amplifier Sensitivity Scale 
Gain Calibration Network 



FIGURE II 



FIGURE 12 
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offset was adjxisted to permit maxinixm flexibility while maintaining a safe 
operating environment for the crystal. A manual offset of 6.8ym was pos- 
sible, insuring that a cavity mode would always be within the biasing 
range. A pair of General Electric Metal Oxide Varistors (V420PA406 Power 
were connected across the output in parallel to the transducer. These 
devices have essentially infinite resistance below a 1000 V input and a 
rapidly decreasing resistance above this value. With these MOV in the circuit 
voltage spikes were prevented from harming the transducer. With these 
adjustments, the transducer was capable of following cavity length fluctua- 
tions for periods of 1/4 - 2 hours without manual adjustment. This period 
is exclusively a function of the room temperature fluctuations. The only 
remaining frequency Instability is due to the oscillation introduced by the 
reference sine wave. This results in a + 4 MHz oscillation about the center 
of the gain curve. 

Cavity alignment procedure was as follows: 

1. A helium-neon laser was positioned so that it passed through 
the mirror hole and was incident on the grating. The laser plasma tube was 
centered around this optic axis. 

2. The zeroth order reflection off the grating was found by 
mantially rotating the grating around the optic axis and locating a spot 
which did not rotate. , 

3. With the reflected orders off the grating all lying in a 
horizontal plane, the 17th order spot was sent back down the optic axis 
and through the hole in the front mirror. Due to mode impurities, not all 
the laser light passes out of cavity. 


- 52 - 



4. The ring of reflected light off the front mirror was then 
positioned along the optic axis by adjusting the front mirror. 

5. The Brewster angle window holders were secxired to the end of 
the glass tube, and slight readjustments of optical mounts were made to 
align cavity exactly. The laser light was Introduced Into a monochomator 

cSi. 

using a f = 4 cm Ge focusing lens. Resolution of the monochromator was 
sufficient; observed slltwldths were five or six times less than the line 
separations. Lasing on a single line was observed for each laser grating 
setting. If more than one line Is lasing, the Intensity of the adjacent 
lines Is less than 1.5% of the major laser line. 

Because a low maintenance system Is desired, pyroelectric detectors 

which operate at ambient temperatures were chosen. These Barnes T->301 

2 

Detectors were 1 mm In area and capable of handling up to 1 W cw. 
Responsivity was 1.5 pa/W for chopping frequencies below 500 Hz. These 
detectors were found to be susceptible to electrical noise and thermal 
fluctuations In the room. Accordingly, they were shielded and enclosed as 
completely as possible. To Insxire that all of the laser beam was Incident 
on the small detector elements, a pair of 1" f/1 Irtran lenses were used 
to focus the light onto the detectors. These lenses also helped Insulate 
the detectors from air currents. The sisals from these detectors are put 
Into a pair of PAR lock*- In amplifiers. With this configuration the meas~ 
ured laser intensity was stable to 0.2% over a few minutes and 0.7% stable 
over the period between manual adjustment of the vortage bias offset. Thus 
the stabilization of the frequency, as expected, results In a stabilization 
of the laser Intensity as well. The dual channel detection scheme was 
Implemented to eliminate errors due to short term intensity fluctuations. 


- 53 - 



The lock-in amplifier's gain may be calibrated for the least sensitive 
scale, but the error in the accuracy of the gain for the more sensitive 
scales may be + 5%. To permit accurate readings of the detector's 
signal for high attenuations (small, signals), a set of multl-tum poten- 
tiometers were Inserted In parallel with the gain trim resistor of each 
lock-in. With proper switching, these resistors permitted accurate cali- 
bration of each sensitivity scale ranging over a factor of 50. This 
circuit Is shown schematically In Figure 12. Calibration was performed 
using a voltage divider network and the internal reference frequency sine 
wave of each lock-in. One sensitivity scale was calibrated against the 
next, and Interscale gain errors were reduced to a negligible magnitude 
In this manner. Any nonllnearltles In response should be due to detector 
limitations and not detection circuit errors. 

A current-buffered analog divider network was constructed to permit 
ratiolng of the two signal channels. An Intronics M 506 Multiplier and 
904" Analog Devices Power Supply were used ^ee Flgurel^. Because the 
Integrated clruclt was designed for multiplication, the division becomes 
unstable for small denominators. Since Interscale switching errors had 
been eliminated, the Input signals into the divider need never be less than 
4V. The output of the divider network is given by 

Y - 10 Z/X (23) 

where Y Is the ratio in volts, Z the signal from the sample channel lock-in, 
and X Is the signal from the reference channel lock-ln. X must be negative 
in sign. The terminals are labeled In the circuit diagram. Calibration 
of the divider network was performed using stable voltage sources 
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Figure 13 : Dual-Beam Divider Network 

The 741C chips act as current buffers for 
analog divider network. and Out Trim 
permit zeroing and calibration. Z is the numer- 
ator, X is the denominator and Y is the output 
ratio. 15 Voltages (+ and -) must be 
supplied. 



■H5 -15 NUMERATOR DENOMINATOR 



FIGURE 13 
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(Heathklt IP 18) and a 4% digit digital umltlmeter (Data Precision Model 
245). After the output trim potentiometer was zeroed for zero Input 
numerator, the trim pots for. the numerator and denominator were adjusted 
to minimize the error for all Inputs between four and ten volts, 

I Z I < I X |. In this manner the error In the divider network could be 
kept to less than 0.7% for all Input signals, and less than 0.2% for 
denominators larger than 8.5 V. The ratio of the laser signal In the 
two beams displayed on the digital multimeter was found to be stable to 
+ 0.8% and reproducible from day to day to greater than 2%. Because 
of the day to day variation, all measurements were made cell In - cell put 

9 

Variable density beam attenuators (polyethylene film) were used to adjust 
for Intensity variation In the different laser lines. These permitted 
the reference channel lock- In to observe a nearly full scale signal at 
each frequency, and as a result, divider circuit errors were minimized. 

The cell used for the absorption measurements was 50 cm long, 

5 cm In diameter, and was fitted with NaCl windows. In this detection 
scheme, all laser output frequencies are colllnear and no variation of 
the ratio of the two channels with frequency was expected or observed. 

The absorption of the empty cell was measured at each wavelength. The 
measurements were found to have a standard deviation of 1.5% (108 values). 
Some attenuation variation with frequency was observed, and remains unex- 
plained. Great care was taken to reposition the absorption cell con- 
sistently. It has since been learned that the pyroelectric detectors 
exhibit a non-linear response acorss the face of the detector, and that 
slight variations in angle and position of the beam incident on the 
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fociislng lens might result la slgaal variations. This change could 
result from thermal expansion of the steel beam supporting the detector or 
slight changes In the spatial position of the laser beam. 


C. Absorption Coefficient Determinations 

Absorption coefficients were determined from equation (21) : 

ONL » eCL « A - -fin (I/I ) , -ln(I/I ) 

( o'cell & sample o cell 

The reproducibility In A values Is listed in Table XIII for various amounts 
of attentiatlon by the sample. It Is seen that large attenuations give rise 


I/Io * 


Tablq XIII 

Data Precision for Various Sample Attenuations 

T3T>ical A with 95% Confidence 
a** Limi ts for Normal Distribution 

(2 a ) 


No. of Points 


0.5-0. 3 

27 

.013 

2.00 ± 

.026 

(1.3%) 

0.15-0.11 

33 

.013 

4.00 ± 

.026 

(0.7%) 

0.07-0.04 

72 

.03 

6.00 ± 

.06 

(1.0%) 

0.025-0.007 

13 

.03 

9.00 ± 

.06 

(0.7%) 


* I/I -(I/I ) , , 

o o sample & cell 


/ 


^^^^o^cell 


** 




to A values which are more reproducible than small attenuations. For small 
at t enuatlons the two terms In parentheses are of conq>arable magnitude and 
the uncertainties In each contribute to the uncertainties In A» while for 
larger attenuations the first term will be dominant and the errors In A 


will be due excliislvely to that term alone. For very small attenuations 
the uncertainties in the two ratios are on the order of their difference 
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and A is conqyletely uncertain. Because the gain of each sensitivity scale 

of the lock~ln was calibrated, little loss of accuracy occurred for small 

signals. Thus an effort was made to obtain large attenuations for as many 

laser wavelengths as possible for each gas» in order to achieve n'g-«-CTmini 

accuracy and precision. The uncertainty in A was estimated to be 3Z for 

l/I between 0.5 and 0.9w 
o 

Reagent grade chemicals were used whenever possible for these deter- 
minations. All samples were further purified by a freeze-pua^thaw 
cycle repeated several times before use. All pressures were measured on 
a dlbutyl phthalate vacuum manometer. Readings were taken with a 
cathetometer and. are accurate to + 1.5% or 0.04 Torr, whichever is 
greater. The manometer response was found to linear and accurate when 
compared with a Hg McLeod gauge. Absorption in the manometer was found 
to be negligible during pressure measurements, and the manometer was 
punq)ed on thoroughly after each use. Pressures up to 35 Torr could be 
measured using the manometer ( p» 1.046 g/cc, 25*^c) . The absorption 
coefficients of the fifteen gases listed in Table XIV were measured 
at the 52 strongest lines of the CO 2 laser. These are the P(6) - P(36) 
and R(6) - R(32) lines of the 00°1 - 10°0 transition and the 
P(8) ■- P(34) and R(10) - R(24) lines of the 00°1 - 02°0 transition. 
Measurements were made both on the low pressure gas sample alone, 
and on the atmospheric pressure mixture with air. A few minutes were 
allowed for mixing of the gas with the air. fk rough calculation shows 
this to be sufficient; see Appendix These coefficients are 

presented graphically in Figures 14-19. Table XV lists the absorption 
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Figure 14 ; Laser line absorption coefficients for 
low pressure trace contaminants ; acetonitrile 
(34.8 Torr); benzene (16.4 Torr); cyclohexane 
(34.7 Torr); 1,2-dichloroe thane (28.9 Torr); 
and ortho-dlchlorobenzene (0.4 Torr). Values 
in meters/Newton ( 1 m/N“ 1 . 33322 Torr ^cm ^) . 

Figure 15 ; Laser line absorption coefficients 
as in Figure 14 for ethyl acetate (5 Torr) ; 
Freon-12 (0.8 Torr); Freon-113 (3 Torr); furan 
(4 Torr); and Isopropanol (10.6 Torr). 

Figure 16 : Laser line absorption coefficients 

as in Figure 14 for methyl chloroform (3.36 Torr); 
methyl ethyl ketone (29.4 Torr); ^-butanol 
(4 Torr); vinyl chloride ((•) 4.1 Torr; (+) 4.1 
TorrS 756 Torr air); and iodopropane (12.3 Torr). 
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.Figure 17 : Laser line absorption coefficients 

for trace contaminants In a mixture with air. 
The partial pressures of contaminant gases 
are the same as In Flgtires 14-16; the total 
pressure Is 1 atm. Acetonitrile; benzene; 
cyclohexane; 1,2-dlchlo roe thane; and ethyl 
acetate. 

Figure 18 ; Laser line absorption coefficients 
as in Figure 17 for Freon-12; Freon- 113; furan; 
and Isopropanol. 

Figure 19 ; Laser line absorption coefficients 
as In Figure 17 for methyl chloroform; methyl 
^ethyl ketone,; t-.butanoi); -yinyl ^chloride 
Torr alone and (+) 4.1 Torr + 756 Torr alr)l 
and lodopropane. 
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coefficients of the gas alone. Table XVI gives the atmospheric mixture 
absorption coefficient values. An uncertainty of + 3% (95% confidence 
limits) exists for all values of absorption coefficients greater than 
4 X 10 ^ m/N. This is due to the combined errors in A and the concentra- 
tlon measurements . It was experimentally determined that fluctuations in 
the value of the ratio give rise to an uncertainty In e of 1 x 10 ^ m/N, 
when A is srull even for large concentrations over 30 Torr. (See equation 
( 21 ).) 

It must be noted that the self-broadening of these gases will present 
a problem, possibly even in the atmospheric mixture. Thus, even though 
the absorption coefficients were determined over a range of pressures for 
each gas in order to obtain good accuracy for both strong and weak extinc- 
tion coefficients, the low pressure measurements were used whenever possible 
High pressure of the sample gas was needed only for frequencies at which 

there was weak attenuation. Self-broadening values have been determined 

221 

for some molecules and are larger than air broadening values. Care must 
be exercised because self broadening could result in overestlmatlon of the 
absorption coefficients and, consequently, underestimation of trace con- 
taminant concentrations. This problem could have been circumvented if a 
longer cell or multiple-pass cell were used for these measurements. A more 
systematic study of the pressure broadening in these gases would be desir- 
able, and these concerns form the basis for the next chapter. For such 

studies, a broadly tunable laser source is preferable to a line-tunable 
222 

CO^ laser. 

It was hoped that use of absorption coefficients taken with low sample 
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Table X7 

Absorption Coefficients of Trace Gases at Low Pressure (m/N)* 


Laser 

ace to- 

benzene 

cycle- 

ethyl 

Freon 12 

Freon 113 

Line 

nitrile 


hexane 

acetate 



R(24) 

.000074 

.000018 

.000008 

.00276 

.00730 

.00068 

R(22) 

.000099 

.000016 * 

.000000 

.00309 

.00690 

.00066 

R(20) 

.000117 

.000019 

.000001 

.00314 

.00470 

.00064 

R(18) 

.000077 

.000017 

.000002 

.00341 

.00353 

.00063 

R(16) 

.000662 

.000009 

.000005 

.00368 

.00347 

.00064 

R(14) 

.000090 

.000010 

.000000 

.00408 

.00261 

.00077 

R(12) 

.000076 

.000006 

.000004 

.00454 

.00222 

.00099 

R(10) 

.000087 

.000021 

.000008 

.00530 

.00161 

.00106 

1 P( 8) 

.000066 

.000431 

.000069 

.0123 

.00023 

. .0107 

1 P(10) 
1 P(12) 

.000125 

.000670 

.000078 

.0120 

.00018 

.0119 

.000126 

.000470 

.000109 

.0119 

.00021 

.0141 

P(14) 

.0^103 

.000458 

.000135. 

.0116 

.00022 

.0159 

o° P(16) 

<N ' 

.000112 

.000520 

: 000125 

.0112 

.00019 

.0177 

7 P(18) 

.000075 

.00118 

.000104 

.0104 

.00019 

.0189 

o-H P(20) 

.000109 

.000359 

.000071 

.00940 

.00018 

.0187 

OP (22) 

.000080 

.000355 

.000095 

.00830 

.00015 

.0174 

P(24) 

.000088 

.000241 

.000146 

.00720 

;.00013' 

.0192 

P(26) 

.000114 

.000308 

.000140 

.00580 

.00009 

.0200 

P(28) 

.000110 

.00126 

.000293 

.00475 

.00020 

.0167 

P(30) 

.000113 

.00217 

.000239 

.00380 

■ .00011 

.0133 

P(32) 

.000082 

.00150 

.000181 

.00296 

.00010 

.0108 

P(34) 

.000070 

.000474 

.000147 

.00227 

.00008 

.0089 

> R(32) 


.000006 

.000006 

.000278 

.00011 

.000468 

R(30) 


.000007 

.000003 

.000257 

.00004 

.000505 

R(28) 


.000006 

.000005 

.000240 

.00005 

.000535 

R(26) 


.000008 

.000000 

.000215 

.00009 

.000530 

* 

1 m/N » 

1.33322 torr“ 

cm ■ 4.1 

X 10 

2 -1 

cm sr at 25 

°C. 
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Table XV 


Absorption Coefficients of Trace Gases at Low Pressures (n/N) 


Laser 

Line 

benzene 

cyclo- 

hexane 

ethyl 

acetate 

Freon 12 

Freon 113 

R(24) 

.000010 

.000002 

.000208 

.00005 

.000570 

R(22) 

.000005 

.000005 

.000187 

.00004 

.000560 

R(20) 

.000007 

.000000 

.000181 

.00004 

.000570 

R(18) 

.000010 

.000001 

.000166 

.00005 

.000550 

R(16) 

.000009 

.000000 

.000158 

.00010 

.000540 

R(14) 

.000006 

.000000 

.000142 

.00005 

.000540 

R(12) 

.000009 

.000000 

- .000140 

.00013 

.000510 

R(10) 

.000005 

.000001 

.000133 

.00007 

.000483 

R( 8) 

.000002 

.000001 

.000128 

.00010 

.000457 

§ R( 6) 
1-1 

.000002 

.000003 

.000120 

.00000 

.000357 

^P( 6) 

.000000 

.000003 

.000142 

.00000 

.000259 

Sp( 8) 

.000002 

.000003 

.000169 

.00010 

.000257 

<-iP(10) 

.000003 

.000002 

.000223 

.00006 

.000257 

o° P(12) 

.000006 

.000005 

.000304 

.00009 

.000277 

° P(14) 

.000005 

.000020 

.000430 

.00013 

.000298 


.000007 

.000011 

.000590 

.00016 

.000326 

°§ P(18) 

.000002 

.000014 

.000700 

.00020 

.000356 

P(20) 

.000011 

.000005 

.000730 

.00033 

.000440 

P(22) 

.000003 

.000006 

.000650 

.00055 

.000580 

P(24) 

.000010 

,000003 

.000660 

.00138 

.000870 

P(26) 

.000006 

.000008 

.000800 

.00365 

.00144 

P(28) 

.000006 

.000009 

.000740 

.0125 

.00224 

P(30) 

.000009 

.000008 

.000840 

.0250 

.00311 

P(32) 

.000007 

.000011 

.000950 

.0220 

.00402 

P(34) 

.000003 

.000017 

.000980 

.0580 

.00560 

P(36) 

.000000 

.000024 

.000960 

.0810 

.00660 
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Table XV 

Absorption Coefficients of Trace Gases at Low Pressure in (m/N) 


Laser 

Line 

furan 

isopro- 

panol 

iodo- 

propane 

tKthyl 

chloro- 

form 

methyl 

ethyl 

ketone 

tert- 

butanol 

vinyl 

chloride 

R(24) 

.000900 

.00275 

.000274 

.00980 

.000463 

.000178 

.00048 

R(22) 

.000800 

.00275 

.000264 

.00820 

.000472 

.000188 

.00043 

R(20) 

.000970 

.00270 

.000239 

.00668 

.000424 

.000177 

.00036 

R(18) 

.000520 

.00279 

.000250 

.00521 

.000431 

.000188 

.00027 

R(16) 

.000725 

.00282 

.000259 

.00422 

.000400 

.000176 

.00034 

R(14) 

.000680 

.00282 

.000242 

.00580 

.000390 

.000181 

.00030 

R(12) 

.00114 

.00298 

.000250 

.00467 

.000377 

.000184 

.00025 

R(10) 

.000840 

.00278 

.000245 

.00476 

.000375 

.000169 

.00023 

g w 

.00119 

.00175 

.000071 

.000870' 

.000270 

.000188 

.00022 

S P(10) 

•H 

.00118 

.00150 

.000068 

.000712 

.000248 

.000185 

^ .00031 

1 

.000796 

.00137 

.000050 

.000546 

.000231 

.000183 

.00028 

H P(14) 

.00122 

.00121 

.000102 

.000449 

.000271 

.000221 

.00031 

^P(16) 

.000774 

.00105 

.000054 

.000296 

.000223 

.000195 

.00090 

P (1'8 ) 

.000699 

.000890 

.000040 

.000222 

'^.000219 

.000216 

^ .00181 

1 P(20) 

.000422 

.000770 

.000055 

.000180 

.000193 

.000222 

.00293 

o^P(22) 

.000920 

,000645? 

.00004# 

.000170^ 

.000205 

.000255' 

.00128 

°P(24) 

.000255 

.000530 

.000069 

.000133 

.000167 

.000266 

.00166 

P(26) 

.000200 

.000442 

.000062 

.000102 

.000171 

.000304 

.00041 

P(28) 

.000182 

.000376 

.000044 

.000099 

.000157 

.000345 

.00145 

P(30) 

.000200 

.000311 

.000088 

.000093 

.000149 

.000388 

.00105 

P(32) 

.000152 

.000272 

.000092 

.000084 

.000121 

.000427 

.00288 

P(34) 

.000157 

.000235 

.000107 

.000070 

.000163 

.000482 

.00041 

R(32) 

.0055 

.000440 

.000040 

.000050 

.000354 

.000125 

.00020 

R(30) 

.0046 

.000515 

.000024 

.000034 ■ 

.000379' 

.000133 

.00029 

R(28) 

*CG157 

.GGu62u 

.000042 

.000026 

.000388 

.000132 

.00119 

R(26) 

.00185 

.000750 

.000025 

.000039 

.000370 

.000105 

.00203 







Table XV 


Absorption Coefficients of Trace Gases at Low Pressure (m/N) 


Laser 

Line 

f uran 

isopro- 

panol 

io do- 
propane 

methyl 

chloro- 

form 

methyl 

ethyl 

ketone 

tert- 

butanol 

vinyl 

chloride 

— R(24) 

.00298 

.000930 

.000028 

.000031 

.000396“ 

.000125 

.00055 

R(22) 

.00108 

.00116 

.000026 

.000039 

.000379 

.000113 

.00051 

R(20) 

.00225 

.00140 

.000022 

.000035 

.000402 

.000110 

.00085 

R(18) 

.00160 

.00165 

.000026 

.000052 

.000385 

.000098 

.00139 

R(16) 

.00090 

.00203 

.000014 

.000050 

.000417 

.000109 

.00051 

R(14) 

.00156 

.00234 

.000007 

.000060 

.000414 

.000094 

.00062 

R(12) 

.00131 

.00263 

.000014 

.000061 

.000418 

.000083 

.00148 

R(10) 

.00126 

.00286 

.000015 

.000056 

.000455 

.000103 

.00110 

R( 8) 

.000920 

.00310 

.000012 

.000051 

.000490 

.000110 

.00053 

g R( 6) 

.000836 

.00316 

.000007 

.000072 

.000515 

.000101 

.00284 

5 P( 6) 

.000074 

.00348 

.000008 

.000080 

.000890 

.000104 

.00116 

§ P( 8) 

.000054 

.00372 

.000009 

.000068 

.000970 

.000128 

.00033 

^ P(10) 

.000034 

.00374 

.000005 

.000056 

.00106 

.000153 

.00053 

P(12) 

.000009 

.00361 

.000016 

.000057 

.00109 

.000169 

.00065 

° P(14) 

.000020 

.00343 

.000011 

.000054 

.00114 

.000225 

.00136 

p(i« 

.000017 

.00316 

.000026 

.000039 

.00118 

.000317 

.00220 

§ P(18) 

.000012 

.00337 

,000015 

.000047 

.00115 

.000450 

.00107 

P(20) 

.000010 

.00308 

.000023 

.000033 

.00121 

.000668 

.00143 

P(22) 

.000006 

.00314 

.000012 

.000038 

.00118 

.001000 

.00814 

P(24) 

.000005 

.00274 

.000006 

.000022 

.00119 

.00137 

.00160 

P(26) 

.000004 

.00267 

.000011 

.000023 

.00114 

.00208 

.00175 

P(28) 

,000010 

.00227 

.000036 

.000030 

.00111 

.00255 

.00102 

P(30) 

.000009 

.00196 

.000045 

.000025 

.00111 

.00305 

.00310 

P(32) 

.000000 

.00167 

.000004 

,000019 

.00105 

.00363 

.00166 

P(34) 

.000003 

.00136 

.000019 

.000028 

.000970 

.00370 

.00499 

P(36) 

.000000 

.00113 

.000017 

.000034 

.000885 

.00347 

.00880 
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Table X7I 

Absorption Coefficients of Trace Gases In an Atm. Mixture with Air (n/N) 


Laser 

Line 

aceto- 

nitrile 

benzene 

cyclo- 

hexane 

1,2-di- 

chloro- 

ethane 

ortho- 

dlchloro- 

benzene 

ethyl 

acetate 

R(24) 

.000073 

.000024 

.000005 

.000004 

.0002 

.00273 

R(22) 

.000075 

.000018 

.000000 

.000011 

.0002 

.00306 

R(20) 

.000100 

.000024 

.000000 

.000003 

.0000 

.00313 

R(18) 

.000075 

.000022 

.000002 

.000012 

.0002 

.00338 

R(16) 

.000412 

.000006 

.000009 

.000005 

.0000 

.00365 

R(14) 

.000084 

.000001 

.000000 

.000009 

.0000 

.00402 

R(12) 

.000081 

.OOOOU 

.000003 

.000008 

.0000 

.00450 

R(10) 

.000092 

.000014 

.000000 

.000004 

.0000 

.00530 

P( 8) 

.000091 

.000438 

.000071 

.000009 

.0004 

.0122 

§ P(10) 

.000122 

.000510 

.000086 

.000013 

.0006 

.0120 

P(12) 

.000128 

.000600 

.000102 

.000015 

.0009 

.0U8 

i P(14) 

.000116 

.000660 

.000113 

.000030 

.0035 

.0117 

^ P(16) 

.000148 

.000680 

.000110 

.000022 

.0037 

.0111 

o® P(18) 

.000097 

.000660 

.000106 

.000019 

.0059 

.0104 

® P(20) 

.000119 

.000580 

.000082 

.000031 

.0026 

.00940 

P (22.) 

.000097 

,.000460 

^^000079 

.000018 

.0076 

.00830 

§ P(24) 

.000114 

.000372 

.000139 

.000041 

.0067 

.00700 

P<26) 

.000127 

.000400 

.000150 

.000035 

.0056 

.00590 

P(28) 

.000096 

.00102 

.000411 

.000034 

.0076 

.00470 

P(30) 

.000U4 

.00210 

.000253 

.000057 

.0074 

.00375 

P(32) 

.000104 

.00111 

.000153 

.000077 

.0064 

.00296 

P<34) 

.000095 

.00069 

.000123 

.000088 

.0054 

.00227 

g R(32) 

.000022 

.000017 

.000008 

.000065 

.0002 

.000277 

®g R(30) 

.000014 

.000009 

.000006 

.000056 

.0000 

.000254 

Ja R(28> 

.000015 

.000008 

.000004 

.000058 

.0002 

.000237 

§ R(26) 

.000012 

.000010 

.000008 

.000054 

.0000 

.000216 





Table XVI 


Absorption. Coefficients of Trace Gases In Atm. Mixture with Air (m/N) 


Laser 

Line 

aceto- 

nitrile 

benzene 

cyclo- 

hexane 

1,2-di- 

chloro- 

ethane 

ortho- 

dichloro- 

benzene 

ethyl 

acetate 

R(24) 

.000019 

.000012 

.000006 

. .000047 

.0000 

.000205 

— R(22) 

.000009 

.000009 

.000007 

.000044 

^.0000 

.000194 

R(20) 

.000016 

.000010 

.000003 

.000048 

.0003 

.000175 

R(18) 

.000002 

.000010 

.00005 

.000039 

.0000 

.000161 

R(16) 

.000009 

.000013 

.000004 

.000045 

.0000 

.000160 

R(14) 

.000004 

.000008 

.000002 

.000050 

.0000 

.000142 

R(12) 

.000000 

.000015 

.000004 

.000050 

.0000 

.000140 

R(10) 

.000008 

.000008 

.000006 

.000064 

.0000 

.000131 

R( 8) 

.000012 

.000007 

.000001 

.000081 

.0005 

.000126 

§ R( 6) 

.000012 

.000000 

.000003 

.000089 

.0004 

.000120 

ap( 6) 

.000001 

.000000 

.000004 

.000197 

.0000 

.000162 

§ P( 8) 

.000009 

.000003 

.000004 

.000215 

.0001 

.000174 

^P(IO) 

.000011 

.000007 

.000008 

.000218 

.0003 

.000217 

P(12) 

.000002 

.000006 

.000013 

.000203 

.0000 

.000302 

S p(14) 

.000013 

.000008 

.000019 

.000175 

.0003 

.000428 

: p(i6) 

.000023 

.000001 

.000013 

.000171 

.0006 

.000580 

§ P(18) 

.000025 

.000008 

.000010 

.000173 

.0000 

.000700 

>^P(20) 

.000045 

.000006 

.000004 

.000200 

.0006 

.000730 

P(22) 

.000048 

.000012 

.000009 

.000208 

.0000 

.000650 

P(24) 

.000058 

.000005 

.000008 

.000204 

.0006 

.000660 

P(26) 

.000085 

.000007 

.000009 

.000197 

.0002 

.000790 

P(28) 

.000090 

.000015 

.000011 

.000187 

.0006 

.000740 

PC30) 

.000098 

.000006 

.000011 

.000163 

.0004 

.000840 

P(32) 

.000116 

.000006 

,000012 

.000129 

.0002 

.000950 

P(34) 

.000099 

.000000 

.000016 

.000118 • - * 

.0000 

.000970 

P(36) 

.000106 

.000005 

.000020 

.000100 

.0002 

.000950 


- 75 - 






Table XVI 


Absorption Coefficients of Trace Gases in an Ata. Mixture with Air (m/N) 


Laser 

Line 

Freon 12 

Freon 113 

furan 

isopro- 

panol 

methyl 

chloro- 

form 

methyl 

ethyl 

ketone 

R(24) 

.00820 

.000680 

.00111 

.00273 

.00980 

.000469 

R(22) 

.00670 

.000680 

.00116 

.00274 

.00830 

.000478 

R(20) 

.00510 

.000640 

.00121 

.00277 

.00680 

.000519 

R(18) 

.00401 

.000640 

.00124 

.00278 

.00531 

.000436 

R(16) 

.00352 

.000640 

.00117 

.00277 

.00434 

.000397 

R(14) 

.00289 

.000780 

.00108 

.00284 - 

.00593 

.000401 

R(12) 

.00209 

.001000 

.00109 

.00297 

.00472 

.000384 

R(10) 

.00165 

.00106 

.000950 

.00284 

.00494 

.000377 

g P( 8) 

.00033 

.0108 

.00121 

.00176 

.000915 

.000277 

2 P(10) 

.00023 

.0120 

.00118 

.00155 

.000630 

.000254 

g P(12) 

.00029 

.0137 

.00111 

.00137 

.000528 

.000226 

^ P(14) 

.00029 

.0162 

.00107 

.00123 

.000441 

.000268 

o pa6) 

.00036 

.0178 

.000880 

.00107 

.000299 

.000223 

g P(18) 

.00020 

.0190 

.000760 

.00092 

.000224 

.000220 

1 P(20) 

.00019 

.0189 

.000562 

.00077 

.000182 

.000199 

®o P<22) 

.00024 

.0175 

.000452 

.00065 

.000168 

.000204 

Z P(24)' 

.00016 

.0195 

.000336 

.00054 

.000138 

.000165 

P(26) 

.00019 

.0202 

.000241 

.000458 

.000104 

.000172 

P(28) 

.00016 

.0169 

.000182 

.000388 

.000116 

'.000161 

P(30) 

.00010 

.0134 

.000153 

.000329 

.000097 

.000148 

P(32) 

.00013 

.0109 

.000134 

.000282 

.000090 

.000120 

P(34) 

.00008 

.0088 

.000129 

.000250 

.000073 

.000155 

^ R(32) 

.00009 

.000468 

.00393 

.000440 

.000050 

.000353 

®o R(30) 

.00009 

.000497 

.00413 

.000520 

.000035 

.000381 

R(28) 

.00002 

.000535 

.00407 

.000630 

.000024 

.000389 

§ R(26) 

.00006 

.000556 

.00369 

.000770 

.000036 

.000368 
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Table XVI 

Absorption Coefficients of Trace Gases in Atm. Mixture with Air (m/H) 


Laser 

Line 

Freon 12 

Freon 113 

furan 

iospro- 

panol 

methyl 

chloro- 

form 

methyl 

ethyl 

ketone 

R(24) 

.00014 

.000563 

.00356 

.000940 

.000031 

.000396 

R(22) 

.00005 

.000573 

.00304 

.00116 

.000039 

.000376 

R(20) 

.00011 

.000569 

.00271 

.00140 

.000040 

.000405 

R(18) 

.00001 

.000556 

.00279 

.00170 

.000053 

.000384 

R(16) 

.00014 

.000531 

.00227 

.00200 

.000057 

.000412 

R(14) 

.00004 

000531 

.00157 

.00232 

.000059 

.000416 

R(12) 

.00004 

.000517 

.00120 

.00265 

.000062 

.000416 

R(10) 

.00013 

.000482 

.000970 

.00292 

.000058 

.000460 

R( 8> 

.00009 

.000451 

.000709 

.00314 

.000053 

.000497 

1 R( 6) 

.00004 

.000412 

.000538 

.00324 

.000073 

.000515 

1 p( 6) 

.00006 

.000259 

.000063 

.00352 

.000070 

.000880 

1 P( 8) 

.00003 

.000256 

.000048 

.00370 

.000067 

.000970 

P(10) 

.00013 

.000258 

.000037 

.00376 

.000055 

.00106 

°o 

.00009 

.000275 

.000025 

.00362 

.000060 

.00110 

7 P(l^) 

.00017 

.000302 

.000020 

.00342 

.000051 

.00115 

o'^ P(16) 

.00014 

.000327 

.000020 

.00319 

.000045 

.00117 

° P(18) 

.00016 

.000364 

.000012 

.00316 

.000043 

.00115 

P(20) 

.00035 

.000450 

.000014 

.00303 

.000038 

.00120 

P(22) 

.00059 

.000597 

.000010 

.00298 

.000034 

.00119 

P(24) 

.00151 

.000910 

.000009 

.00280 

.000024 

.00120 

P(26) 

.00367 

.00146 

.000008 

.00261 

.000027 

.00115 

P(28) 

.00870 

.00232 

.000011 

.00232 

.000029 

.00112 

P(30) 

.0180 

.00321 

.000007 

.00198 

.000023 

.00111 

P(32) 

.0308 

.00407 

.000005 

.00167 _ 

.000028 

.00106 

P(34) 

.0429 

.00560 

.000005 

.00138 

.000031 

.000970 

P(36) 

.0463 

.00660 

.000004 

.00114 

.000030 

.000890 
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Table XVI 

Absorption Coefficients of Trace Gases in Atm. Miztxire with Air (m/N) 


Laser 

tertiary 

vinyl 

Laser 

tertiary 

vinyl 

Line 

butanol 

chloride 

Line 

butanol 

chloride 

R(24) 

.000170 

.00000 

R(24) 

.000130 

.00154 

R(22) 

.000202 

.00022 

R(22) 

.000120 

.00093 

R(20) 

.000170 

.00000 

R(20) 

.000117 

.00139 

R(18) 

.000195 

.00015 

R(18) 

.000100 

.00104 

R(16) 

.000171 

.00004 

R(16) 

.000106 

.00105 

R(14) 

.000187 

.00000 

R(14) 

.000100 

.00126 

R(12) 

.000186 

.00000 

R(12) 

.000086 

.00096 

R(10) 

.000163 

.00002 

R(10) 

.000100 

.00175 

e P( 8) 

.000189 

.00017 

R( 8) 
c 

.000108 

.00099 

H P(10) 

.000189 

.00023 

5 R( 6) 

U 

.000094 

.00247 

•H . 






a P(12) 

.000189 

.0003L 

s « 

.000101 

.00147 

H P(14) 

.000227 

.00046 

H P( 8) 

.000130 

.00114 

"" P(16) 

.000200 

.00069 

P(10) 

.000162 

.00130 




o 



P(18) 

.000220 

.00085 

°o P(12) 

.000163 

.00136 

o 



iH 



j P(20) 

.000229 

.00105 

1 P(14) 

.000223 

.00128 

P(22) 

.000257 

.00140 

o'^ P(16) 

.000322 

.00149 

o 



o 



° P<24) 

.000263 

.00128 

® P(18) 

.000438 

.00172 

P(26) 

.000306 

.00140 

P(20) 

.000666 

.00224 

P(28) 

.000351 

.00145 

P(22) 

.000990 

.00884 

P(30) 

.000384 

.00133 

P(24) 

.00144 

.00321 

P(32) 

.000432 

.00123 

P(26) 

.00207 

.00233 

P(34) 

.000480 

.00098 

P(28) 

.00255 

.00322 

.00180 

.00193 




Pt3u; 

^ R(32) 
o 

.000110 

.00068 

P(32) 

.00364 

.00169 

®o R(30) 

pH 

.000137 

.00077 

P(34) 

.00380 

.00169 

^ R(28) 
o 

.000135 

.00107 

P(36) 

.00353 

.00317 

2 R(26) 

.000113 

Ann TO 
• %/%/w / w 
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pressures would e^xlbit self-broadening which was negligible coiq>ared 
with atmospheric broadening, and accurate detection of each gas In a 
mixture would be possible. This Is discussed further later In the 
chapter. 

From these absorption coefficients, the sensitivity of a laser 
absorption method for detecting each trace gas contaminant In the 
atmosphere could be estimated. The minimum detectable concentration 
(MDC) of each gas was set at the level for which there Is a 20% uncer- 
tainty In A at the gases' most strongly attenuating frequency. These 
MDC values are listed In Table XIV. It may be seen from the Table that 
single gases may be detected at the 1-200 ppm level for a 10m path. 

The gases to be used In preparing a mixture were purified as above. 
After an amount of each gas reached equilibrium in the gas manifold 
and Its pressure was measured, the cell was closed off from the manifold 
and the gas In the cell was cryopusq>ed Into the sldearm bulb of the 
cell, using liquid nitrogen. This bulb was of known volume relative to 
the volume of the absorption cell. Care was taken to Insure all the 
contaminant was condensed Into the bulb. After all the gases were 
stored In the bulb In this manner, the bulb was warmed to room tempera- 
ture and the gases were allowed to expand into the absorption cell and 
mix for several hours. Air was then added to the cell to bring the mix- 
tures to atmospheric pressure. Calculations of the times necessary for 
diffusion are given in Appendix II. " * 

Mixtures were prepared In such a way as to find the effect of the 
following conditions on the accuracy of the concentration determination 
procedure: 
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1. The variation of the absorption for similar mlxtures;l.e. , 
mlxtxires with the same relative concentrations of the gases, but different 
total pressures of the contaminant mixture; 

2. The variation of the relative concentrations of two gases 
with very similar absorption spectra; 

3. The presence of a relatively high concentration of a broadly 
absorbing gas. 

D. Mixture Analysis — «« • 

223 

Several different approaches to the solution process are possible. 

A best fitting of the calculated concentration by a linear vector space 

132 

analysis for a single species has been performed. Here four wavelengths 
of a CO 2 laser were chosen to permit maximum discrimination for 0^ detec- 
tion. A least-squares fitting of a high resolution continuous spectrum 

91 

has been done for up to 25 contaminant species. This method relies on 
determining the concentration of each species from a distinguishing feature 
or set of features In the absorption spectrum. The method chosen here 
to determine the concentration of each gas present In the prepared mix- 
tures involved a least-squares fitting of the absorption data prior to 
solution. As a result, the magnitude of the errors due to experimental 
fluctuations was minimized, and greater stability during matrix manlptila- 
tlon was achieved. 

Because only as many slmultaneoiis equations were required as 
unknowns, fifty-two absorption measurements on a twelve component mixture 
form an over determined system. Early atteoq>ts at choosing appropriate 
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132 

wavelengths for discrimination as In for even a three component system 

resulted in low accuracy and a large dependence of the solution on the 

choice of wavelengths. It would seem that the linear vector space 

132 

analysis used by Craig worked sufficiently well for the simple one- 
component system, but that poor stability would result for a multl- 
conq>onent mixture using their Iterative approach. As a restilt, a method 
utilizing al T the experimental Information was Implemented. The expected 
absorption at each wavelength. If the absorption due to each species 
“ pfSSent is additive, would be given by 






,X 

,th 


(24) 
the 


where y^ Is the expected absorption at the 1 wavelength, Is 

experimentally determined absorption coefficient of gas j at wavelength 1, 
and Cj Is the "unknown" concentration of gas j. The experimentally measured 
absorption values of the mixture, A^, are subject to the uncertainties and 
fluctuation errors'dlscussed above, and as a result will generally be 
different from y^ values. The best set of solutions for tbe concentrations 
of the contaminants occurs when the differences between and y^ values 
are smallest (but not necessarily zero) . This solution Is achieved when 
the Slim of the squares of the residuals between A^ and y^ are minimized with 
respect to variation of the solutions c^ , giving 


^ 

— » 0 ( j = 1,...,N) ; n “ I (y^. - A )‘ 

3Cj 1-1 ^ ^ 


(25) 


For this system N « 12, X 52. Substituting equation (24) Into equations 
(25) and performing the partial differentiation results In a set of N 
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simultaneous equations after rearrangement. 


^j2 ^j3 V 


( 26 ) 


j - 1 , 2 , ... . N . 


where 


Jk 


X 

I e 

1=1 


Ij ^Ik 


and 


I K / 2 > 
1=1 


J .i, ‘■ij "i 


(27) 


It Is seen that each and dj matrix element Is the sum of fifty-two 

products. In this manner all the Information Is collapsed Into an 

N-dlmenslonal matrix. In matrix notation equation (26) may be written 

EG “ Q where E = g g and D = £ A /I . Because the F matrix Is 

symmetric positive definite, elimination methods requiring non-vanishing 

positive pivots may be employed. In particular, the triangular decomposl- 

224 

tion method of Cholesky which breaks F into two trlangtilar matrices 
which are transpositions of one another, introduces very little mathe- 
matical Instability Into the solution process throxigh minimization of 
the rounding errors In the matrix elements calculated during the 


decomp.qsltiou. The ’-solutions : are '•then ;^found';by 'back-substitution. -A 

fuller explanation and an example are given in Appendix III. Methods 

requiring matrix Inversion for solution are far more unstable. They 

are partlctilarly imstable for near-vanishing pivot elements. This would 

be the case if the sum of the absorption coefficients for two gases were 

X 

I 

1=1 

matrix Inversion methods were unable to reach stable solutions for the 
twelve component mixtures in these experiments. The coiiputer solution 
program, DPELZC, Is given In Appendix IV. An IBM System/360 Model 67 


£ ) ^ 
^ik^ 


)• 


We noticed that. In fact. 


similar 


r ^ 

I ( 

1 1-1 




2 - 
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co^uter was tised in this analysis on a time-sharing basis. The memory 
storage of the PDP 8/L mini-computer In the laboratory was Insxifflclent 
for these large matrices, md the word length did not allow sufficient 
precision during computation. 

Occasionally gases were omitted from the mixture and If the cal- 
culated concentration of the gas was negative or zero, that species was 
removed from the calculation process, redticlng the nuoiber of slmultaneotia 
equations by one. If the calculated concentration of an omitted gas 
was non-zero, or If the concaitratlon of a trace gas present In the mix- 
ture was calculated to be zero, the error In the solution was taken as a 
total uncertainty In the presence of the contaminant at that level. Gen- 
erally, as the concentration of any gas decreased, the percent error 
between the calculated and known Input values Increased even for similar 
mixtures. This Increase Is expected because the relative contribution of 
that gas to the absorption by the mixture Is decreased and closer to the 
noise level of the detection system. However, these errors were a function 
of the Interference effects In each particular mixture and not just the 
contaminant's concentration. For example, benzene and cyclohexane have 
very similar absorption coefficient distributions, as listed In Table XVII. 
Benzene absorbs more strongly In this region and Its maxlnnm absorption 
Is shifted by one CO 2 laser line toward lower frequency. The MDG of 
cyclohexane Is Increased when benzene Is present In an amount sufficient 
to dominate the relative absorption. Under these conditions, the benzene 
MDC Is decreased below the value expected from concentration considerations 
alone. The concentration effect Is predominant over most of the relative 
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Table XVII 


The Absorption Coefficients in m/N of Benzene and 
Cyclohexane in the P Branch of the (00°1 - 02°0) 
Transition * 


Transition 

Frequency (cm ^) 

(Benzene) 

(Cyclohexane) 

P (8) 

1057.3002 

.00043 (1) 

.000070 (2) 

P (10) 

1055.6251 

.00050 (1) 

.000084 (2) 

P (12) 

1053.9235 

.00059 (1) 

.OOOiqgsv^):^ . 

P (14) 

1052.1956 

.00065 (1) 

.000111 (^T 

P (16) 

1050.4413 

.00066 (1) 

.000108 (3) 

P (18) 

1048.6608 

.00065 (1) 

.000104 (3) 

P (20) 

1046.8542 

.00057 (1) 

.000080 (3) 

P (22) 

1045.0217 

.00045 (1) 

.000077 (3) 

P (24) 

1043.1633 

.00036 (1) 

.000136 (4) 

P (26) 

1041.2791 

.00039 (1) 

.000147 (4) 

P (28) 

1039.3693 

.00100 (1) 

.000403 (8) 

P (30) 

1037.4341 

.00206 (2) 

.000248 (5) 

P (32) 

1035.4736 

.00108 (1) 

.000150 (4) 

P (34) 

1033.4881 

.00068 (1) 

.000121 (4) 

* The absorption coefficients 

at all other 

laser frequencies are 

near zero 

for both gases. 





concentration range investigated, and th\is this method of solution is suffi- 
ciently sensitive to easily distinguish the two gases. 

The concentrations of the gases at which the difference between cal- 
culated and known values becomes sizable for the twelve conq>onent mixtures 
are given in Table XV. It is seen that the presence of other absorbers 
tends to. Increase the MDC value of each contaminant, general to about 
five times that of the isolated gas. Mixtures containing only a few coot- 
ponents should have MDC values intermediate to those of the twelve com- 
ponenW^^lpre and single gas MDC levels. The errors listed in Table XIV 
are the observed differences for the gas at that concentration in a mix- 
ture. Interference effect do not permit the MDC levels to be set as 
uniformly as for the isolated gases. 

Another interference effect in the mixture may be used to increase 
the sensitivity of detection. When all the gases are present in low con- 
centrations, the absorption of the laser light by the mixture is small, 
and the values are largely uncertain due to baseline fluctuations. 

The addition of a broadly absorbing gas decreases the effect of the 
fluctuations and actually increases the sensitivity of detection of the 
other trace species. For example, isopropanol exhibits two broad absorp- 
tions in the region of the spectrum covered by the laser ^ee Figure 1 ^. 

Two mixtures of twelve gases were prepared containing concentrations of 
the gases at or below their MDC values determined above. One contained 
an amount of isopropanol near its MDC value, and the bother contained 
an amount roughly nine times MDC, resulting in 80-95% absorption of the 
light at many frequencies. The calculated concentrations are conq>ared in 
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Figure 20 ; The absorption coefficients of 
Isopropanol at the CO 2 laser frequencies. 

The frequencies at which the following gases 
have strong absorption coefficients are 
indicated: acetonitrile (A), 1, 2-dichloro- 

'ethane (D), ethyl acetate (E), Freon 12 (12) 
Freon 113 (113) , and vinyl chloride (V) . 






Table XVIII. A large increase in sensitivity resulted for those contami- 
nants which have their strongest absorptions in a region where the addition 
of Isopropanol greatly increases the total absorption. Acetonitrile, 

1,2 dichloroethane, and vinyl chloride are in this category ^^e Figure 2 ^. 
The frequencies at which ethyl acetate. Freon 12, and Freon 113 have their 
strongest absorption coefficients are also indicated in the Figure. It 
is seen that the total absorption at those wavelengths wotild not be 
greatly Increased by the addition of isopropanbl, and little Improyement 
is observed in Table XVIII. . 

Normal atmospheric attenuation over a 10 m path due to water vapor 

131 133 

and carbon dioxide is less than 0.2% for all frequencies. ’In 

225 

closed environments, both these gases may acctunulate to high levels. 

If their absorption is significant, they may be treated as additional 

atmospheric contaminants and an attenpt could be made to quantitatively 

detect them. Water vapor would give rise to, at most, 0.3% absorption 
132 

for a 10 m path, and hence is imdetectable since experimental measure- 
ment fluctiiations are on the order of ' 1%. Carbon dioxide, on the other." 
hand, might be present in sufficient quantities to restilt in 10% absorption 
at some laser frequencies, and thus could be quantitatively detected. In 
these closed environments, ambient CO^ may serve as a broad absorber and 
enhance sensitivity. The additional errors introduced into the 'solution 
process would be small and- the accuracy of the detection method would then 
be independent of cabin atmosphere composition. The presence of unexpected 
trace absorbers will greatly perturb the ealcvilated solutions. As a 
result, all species which might occur in the closed environment must be 
enumerated and the appropriate absorption coefficients stored in memory. 
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Table XVIII 


The Effect of a Broadly Absorbing Gas. on the MDC 
Values in ppm * 


Ratio of the calculated 

concentration to the MDC of gas 


Gas 

known value for the 
mixture with: 

All Gases Addition of 
near MDC Isopropanol 

MDC of gas 
from 

Table XV 
1 = 10 m 

due to 
presence of 
Isopropanol 
1 =a 10 m 

acetonitrile 

2.62 

1.34 


400 

(33%) 

140 

(33%) 

benzene 

1.82 

1.10 


200 

( 6%) 

50 

(10%) 

cyclohexane 

1.54 

0.64 


167 

(50%) 

130 

(30%) 

1 , 2-dichloroe thane 

1.02 

0.52 


1000 

(30%) 

125 

(50%) 

ethyl 








acetate 

2.45 

1.18 


9 

(23%) 

9 

(15%) 

Freon 12 

0.66 

0.85 


5 

(38%) 

9 

(15%) 

Freon 113 

0.89 

1.16 


9.2 

( 7%) 

9 

(22%) 

furan 

1.12 

0.93 


33 

(20%) 

40 

( 7%) 

Isopropanol 

2.64 

1.04 


40 

(10%) 

- 


methyl chloroform 

0.82 

0.95 


9 

(14%) 

20 

( 5%) 

t-butanol 

2.98 

1.65 


73 

(25%) 

47 

(58%) 

vinyl chloride 

0.49 

0.81 


17 

(50%) 

17 

(18%) 

— 

* 

1 N/m^ - 9. 

878 

ppm 





- 89 - 




Elimination of species not actually present will permit eventual matrix 
slnq>llflcatlon. 

The enhanced sensitivity brought about by the addition of a broadly 
absorbing gas verifies the fact that more strongly attenuated laser 
frequencies give absorptions which are more reproducible. Thus one 
might consider weighting the absorptions differentially to try to inq)rove 
the accuracy of the solution method. In this case equation (26) becomes 

e^uieC»e^(DA/il (28) 

where oj is a X x X diagonal weighting matrix. Heavier weighting of 
strong absoirptlons was tried initially, but a decrease in accuracy 
resulted. The least-squares fitting part of the solution process 
tends to treat large absorptions preferentially, and no improvement is 
gained by weighting. More heavily weighting the small or medium absorp- 
tions was also used in an attempt to determine whether or not the least- 
squares fitting process stressed the large absorptions too strongly. 

These solutions were not significantly lnq>roved over the unweighted 
scheme. Apparently the imblased least-squares fit of the data produces 
the best set of solutions and further attesqjts at weighting offer no sig- 
nificant improvement. A weighted computer solution program, W2DPAZC, is 
given in Appendix IV. 

Because the MDC depends on Interferences in the mixtures, a com- 
putational technique was developed which assigns a relative uncertainty 
to the calculated concentration values for each specific mixture. No addi= 
tlonal eiqperlmental data was needed and the extra computer time required 
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is Inslgxilflcant. Such a self-evaluating procedure Is desirable when Inter- 
ference with other absorbers may cause a trace species to become partially 
obscured or much too prominent. 

The first approach to this self-evaltiatlon Involved calculation of 

227 

the residuals between the calculated and exact solutions of the matrix. 

= The Cj solutions were used to calculate y^ value tislag»equatlon (24). 

The difference between y^ and Is taken and this value Is used as the 
value In equation (26) and (27). The D matrix calculated Is the 
residual due to computer rounding and experimental fluctxiatlon errors. 

If the experimental A^ values were exact, then the residual would reflect 
the mathematical stability. The c^ values calculated using the residual 
D matrix are the errors In the concentration values previously calculated. 
The relative Importance of the two sources of error can be estimated as 
follows. If the number of certain digits In the experimental values Is 
q and the word length of the computer is p. It follows that If the uncer- 
talnty in the (p + 1) digit gives rise to a residual then the uncer- 

talnty In the (q + 1) digit gives rise to a minimum error In the concen- 
tration solutions of the order of magnitude 10^^ Thus the error 

due to physical tmcertalntles > 10^^ and the residuals calculated 

by this process, R, are R < (1 + 10 For the IBM System/ 360 

Model 67 computer used In this imalysls, p 7. The standard deviation 
In the experimental absorptions Is less than 0.5% for all reasonable 
absorptions Js&e Table XII^. The D vector is composed of the product of 
these values and the standard deviation In the mean of each dj value Is 


related to that of the A^ values which compose It, 


228 


/T 


0.139 a. = 0.07% 
A 


(29) 
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The conq>uter rounding errors may be considered to be completely random. 

4 

Hence q = 3 and *10 > R. The mathematical errors Introduced by 

the solution process are negligible compared to the errors Introduced 
by the experimental uncertainties. The residuals calculated using this 
technique give rise to an estimate of the uncertainties In the concentra- . 
tlons of 0.03% or greater, d\ae exclusively to the experimental measurement 
fluctuations. 

However, the observed errors are much larger than this and so a finer 
estimate of the effect of the physical errors on the stability and accuracy 
of the solutions Is needed. A method allowing modification of the experi- 
mentally determined absorption values was implemented. After the concen- 
trations are found, these values are used to calculate y^, the expected 
absorption at the first wavelength. This value Is compared to A^ and Aj^ 

Is allowed to vary by up to 2% In the direction of y^^. This new value 
(Aj^’), along with the other 51 unchanged A^ values. Is used to recalculate 
the concentrations, which are. In tiu:n, used to adjust A^. This Is 
repeated for all the wavelengths, and the entire process may be repeated 
xmtll the differences between y^ and the modified A^ are arbitrarily small. 
Typically, fewer than five Iterations were needed to make y^ - A^' < 0.001. 
The mathematical round errors accumulate, however, and are as much as 260 
times greater than for the simple solution. This Is still less than the 
mlnlimnn physical uncertainty expressed above. The final set of concentra- 
tions represent the exact solutions of the equations using modified absorp- 
tions; This minimization of absorption method produced solutions whose 
values were different from the original calculated solutions by only a 
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small amount for gases far above their MDC; i.e., for gases responsible 
-for significant absorption of the light. For gases at concentrations near 
their MDC, the differences were larger by up to a factor of 10 than the 
error in the original calculated solutions. These new solutions are not 
necessarily in the direction of decreased error co9q>ared to the known 
input concentrations and are not meant to improve the accxiracy of the cal- 
culated solutions, but are Intended to provide a better estimate of the 
magnitude of the errors due to the experimental uncertainties for each 
gas in each particular mixture. This method provides a better estimate 
of the uncertainties in the solutions than those calculated from the 
residuals. which form a lower bound for these errors. Observation of the 
minimized absorptions reveals that the strong absorptions are seldom 
changed by as much as 2% after 5 iterations ' (max. 10% change possible) , 
whereas the weak absorptions (I/I^ > 0.4) are greatly altered, with the 
weakest changed by nearly 10%. Thus compounds whose strongest absorption 
coefficients are in regions of low absorption of the mixture may have their 
concentrations drastically altered by this minimization process. This 
effect is noticeable in the solutions for concentrations of compounds that 
are weak absorbers , which are made very unstable by this iteration method. 
The preferential treatment received by strong absorbers in the simple 
solution process, is amplified by this method, which highlights the 
instabilities. 

Finally, it must be noted that no systematic deviation of the calcu- 
lated concentrations from the known values occurred at very low con- 
taminant levels. If self-broadening were significant in the measured 
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absorption coefficients a consistently low concentration value would have 
been calculated. The least-squares fitting procedure emphasizes the 
role of the stronger absorption coefficients which were determined by 
low trace gas concentrations. Thus, it is believed that there is some 
self -broadening contribution to the observed absorption coefficients, 
but the solution procedure used, results in a minimization of these 
effects to the point where no experimental perturbations are observable. 
Nonetheless, a more systematic study of pressure-broadening effects in 
many of these systems would be desirable. For such studies, a broadly 
tunable laser source is preferable to a line-tunable system such as C02> 
The use of such a source permits one to map out a complete profile,* 
rather than just sample a single point on that profile. 


E. Conclusions 

It has been shown that the CO2 laser absorption technique is capable 
of sensitive, real-time detection of trace contaminants. The method per- 
mits -simultaneous detection of ;f if teen or more species at ppm .'levels . 


This sensiti-vlty is largely due to the stability and reproducibility of 
data gathered using the frequency stabilized CO^ laser and dual beam 
detection system. Absorption Interferences have been found not to seri- 
ously degrade performance from the single gas levels. In addition, the 
presence of absorbers , like CO2 > in the ^UIlblent atmosphere may actually 
Increase the sensitivity of detection. The development of a self- 
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varying ambient conditions. 
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The standards set by the Cooanittee on Toxicology of the NAS Panel 
20 

on Air Quality for both short- and long-term expostures are listed in 
Table XV. The MDC's we are able to obtain for five of the gases listed 
In the Table are below the long-term exposure level, and an additional 
four are below the short-term level- In the case of vinyl chloride, the 
’exposure levels have been set particularly low because of the carcinogenic 
properties of this gas. Extending the sampling path from 10 to lOO a would 
permit detection at or below this level. Such a configuration would be 
appropriate for an industrial monitoring situation, in which vinyl chloride 
contamination is most likely to be a problem. We have also seen that addi- 
tion of a broadly absorbing gas can enhance sensitivity to low concentra- 
tions. This might be practicable for samples withdrawn from the atmosphere 
into an enclosed multiple-pass monitoring cell. 

Currently precision is limited by the experimental variations in the 
absorption measurements. Either reduction in the magnitude of these fluc- 
tuations or utilization of a greater number of experimental wavelengths 
would improve calculated concentration precision through the least-squares 
fitting procedure of the absorption data. It is seen from equation (29) 
that the expected deviation in the D matrix is proportional to the inverse 
root of the number of different wavelengths used. For 52 wavelengths and 
15 gases, the mathematical contribution to the error in the calculated con- 
centration is negligible. Increasing the number of wavelengths should 
.permit a corresponding increase in the number of simultaneously detectable 
gases at a reasonable level of precision and sensitivity. One way to accom- 
plish this would be to use several isotoplcally pure CO^ lasers. In this 


; 
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manner several hundred lines spanning the 8.9 ym - 11.4 ym wavelength 
region could be employed. Such a system would be bulky and have large 
power requirements; separate laser cavities and support electronics would 
be needed for each laser. Use of a mixture of all the CO 2 Isotopes Is 
not possible due to quenching. 

Another way in which this could be achieved would be to use continx*- 

o . o . * . 

ously tunable, rather than fixed- frequency (line-tunable) monitoring 

sources . If a given pollutant has sharp spectral features , then the laser 

can be tuned to one of these sharp features, free from interference, and 

good analyses can be obtained. This has been accomplished, for example, 

for NO using a spin-flip Raman laser^® and for ethylene, and 

221 

SO 2 using a tunable diode laser. In spectral regions where all the 
contaminants absorb, the data reduction methods we have employed may be 
used to reduce these absorption measurements to derive concentrations 
without running Into computer limitations. If semiconductor diodes were 
used as the tunable sources, an array of these tiny lasers could be 
powered by the same sx^port eqxilpii»nt, and they could be chosen to 
Insure that all contaminants of Interest would have significant absorption 
at some laser wavelength. 



APPENDIX 


Concentration Determination Procedure 


A. Least-Squares Data Fitting 

Substituting in for y. in equation (25) using equation (24): 


X X f N 

n = I ( “ A )2 = I { ^ e C 2} - A 

i=l i=l'‘ 3=1 ^ J- 

As an example, let X = 4 and N = 3, then 


(Al) 


n - I 

1=1 


(Ell “=1 ^12 “=2 * =^13 “=3 - ^1 


(A2) 


= I ■‘■^12 ^^2 ^ +£^3032)^+ aJ - 2 A^(e^^c^2 +£^2*^2^ 

ei3Cjl)) 

Performing one of the partial differentiations. 


or 


• 5 ^ - 0 = ^ ( 2e ^ c, 2^+ 2e,, £.. c„2^+ 2e .e c,2^- 2A e.,2) 

3c, il 1 il 12 2 13 il 3 i il 

1 1=1 

X X 

^il^il^l ^11^12^2 ^il^i3‘^3^ “ ^il^l^^ (A3) 

l*i i=x 


and using equation (27) 


( fj_j_ Cj_ + fj_2 c^ + fj_3 C3 ) - dj_ 


(A4) 


Repeating the differentiation for each c^ gives rise to an ( N x N ) F 
matrix of the form given in equation (26) . 


Matrix Decomposition 

Normal matrix decomposition products are of the form £ = L U 
where L and U are triangular matrices below and above the diagonal. 
Both L and U matrices need to be stored. Either L or U will not be 
a unit triangle. If the matrix F is symmetric and positive definite. 
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as must be the case for matrices obtained from a least-squares fitting, 
the decomposition may take the form F = L L' where the triangles are 
the transpose of each other. I, is now not a unit triangle. This is 
seen to be equivalent to the deconroosition F = L T L’ , where L is 


a unit triangle and T is diagonal non-unltary, 
where is the square root of T 
be recorded. For this system 


Ik 

if we take L = L T , 
~ -u ~ ’ 

1 . Only L or L' need 


‘ll 

^12 

^IN 

^21 

• ' 

^22 • * * 
• 

^2N 

* 

• 

‘ ^N1 

• 

• 

^N2 

^NN 


r \ 


'N 






(A5) 


The L' triangle is found by the following procedure: 

1) =■ ( ^ 11 * Pi *'11 first row 

for all s > 1. 

2) let » - f^^/ then set f^^ = q^^ + f^^ and 

dj " qjj^ p^ j “ 2, . . . ,N and s » 1, . . . ,N. 

3) 4^2 “ < *22>'‘- “2s ■ *2s' “ 22 ' ”2 ' ^ 2 ' hi 

for all s > 2. • 

4) let - - 'j2'' “ 22 ’ “js " ’j2 “ 2 s “je 

dj « q ^2 ?2 j * 5, ... ,N and s = 2, ... ,N. 

5) continue until triangularlzed. 

Now concentration values may be found by backsubstitution in equation 


V 


c 


P 


As an example consider the system of equations: 
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2 

1 

1 


1 

2 

1 


1 

1 

1 



C-I 


’ 5 ' 


1 

c_ 

a 

6 


2 

c_ 


4 

J 

1 3j 




Solving, using the steps Indicated above. 


1) 


1/y^ 

1/vT 

5/»^ ■ 

2) 

yT 

1/v^ 

H'/l 

• « 

5/71 


1 . 

2 

1 

6 


0 

3/2 

“l72 

7/2 


1 

k ■ 

1 

1 

4 


0 

r 

1/2 

1/2 

3/2 


3) 


5) 


1/y^ 1/v^ 

0 1/572 

0 1/2 1/2 

l/v^ 1/t^ 

0 7572 

0 0 vT73 


5/*^ 

VlJTTl 

3/2 

= L’ 


5/71 

ViJipi 

Kfl 


4) 


1/71 l/y^ 

5!/Z 


0 

7571 h-Zin 

7/27573 


0 

0 1/3 

1/3 . 


By back-substitution * 1 » ^2 ** 
programs are given In Appendix IV. 


The computer solution 


(C . Calculation of Residuals 

Once the solutions have been found, the magnitude of the mathemat 
Iced errors In the solution process may be estimated by the following 
procedure; 

1) Compute = D - F c using double precision. 

2) Solve F 6c = R^ using the L' matrix already calculated. 

3) Set c' " c + 6c, rounding to single length. 

4) Repeat using c' Instead of c In step 1). 
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APPENDIX : COMPUTER CONCENTRATION DETERMINATION PROGRAMS 
Elimination Method with Zeroing of Concentrations 

C \s.oarce foe Rouble _Precision _^imination with Zeroinq of 
C neqative Conceatrations. 

c 

C SOURCE. DPELZe 

C . 

C This proqcam solves a system of LAMBDA eqns. in NGASES 
C unknowns, where LAMBDA>NGASES. LAMBDA may not exceed 
C 55; NGASES may not exceed 15. Both are read as input. 

C Only non-zero inputs need be entered, one element at a 
C time preceeded by its row and column numbers. A row 
C number of 99 siqnals end of data. The proqcam requires 
C the followinq parameters prior to the enterinq of data; 

C LAMBDA: number of wavelenqths for which there is data 
C NGASES; number of qasas in the mixture for this run 
C BIGSTA; the maximum absorption to oe permitted 
C BIGSTE: the larqest absorption coefficient to be permitted 
C All input is checked for validity. 

c 

DIMENSION ABSOfiP (55) ,ABSCOF(55,15) ,C0NC(15) ,ACeSD(15,15) , 
♦ D(15) ,INDEX(15) 

DOUBLE PRECISION ACRED, D,SD M, TEMP, Q 
LOGICAJ. OK 

First clear arrays of data to be read in. Create index. 

DO 20 1=1,55 

ABSOfiP (I) =0.0 
DO 10 J= 1, 15 

ABSCOFd, J) =0,0 
10 CONTINUE 

20 CONTINUE 

DO 25 J=l,15 
CONC(J)=0.0 
IHD£X(J)=J 
25 CONTINUE 

C 
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C Bead control parameters now. — ^ 

BE&D{5,5C0) LAiiBDAyMGASES.BIGSIA. BIGSTf; 
NGASP1=HGASES*1 

IF (.HOT (LAMBDA. GT. 55) .OH. (LAMBDA. LT. 1) ) ) GO TO 30 
IF(. MOT. ( (NGASfiS. GT. 15) .OH. (NGASSS.lt. 1» ) ) GO TO 3d 
IF(.NOT. ( (BIGSTE. LT. 0.0) .OB. ( BIG SIA . LT. 3 . 0) ) ) GO TO 30 
H3ITE(6,510) 

STOP 

If input parameters ware bad, 1ob was aborted. Read in 
arrays ABSOHP first. Set counter of bad data to 0. 

Check data validity before storinq. 

30 HEHfiOR=0 

LIMIT=LAMBDA*1 
DO 35 KAHDNO=1, LIMIT 
H£AD(4,515) I, J, TEMP 
I? (I.EQ.99) GO TO 36 
OK^.TRUE. 

IF( (I.LI . 1) .OR. ( I. GT. LAMBDA) . OB . ( J . NE . 1) ) OK= .FALSE. 
IF (DABS (TEMP) .GT .BIGSTA) OK=. FALSE. 

IF (OK) ABSOHP (I) =TEMP 
IF(.HOT.OK) WRITE (6,530) I, J, TEMP 
IF(.NOT.OK) N£fiHOR=NERHOR+1 

35 CONTINUE 

If DO is satisfied, there were too many absorptions read. 

WRITE (6,540) KARDNO 
STOP 

Head in absorption coefficients. NEEHOE is cumulative. 

36 LIMIT=NSASSS*(LAMBDA) +1 
DO 39 KARDON= 1, LIMIT 

READ (7 ,520) I, J, TEMP 
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IF (I.EQ.99) GO TO 4C 
OK=.TBUE. 

IFCd.LT. 1) .OB. (I.GT.LAaBDA))OK=.FALSE. 

IF ( (J. LT. 1) .OR. (J.GT. NGASES) ) 0K=. FALSE. 
IP(DABS(TEflP) . GT. BIGSTE) 0K=. FALSE. 

IF (OK) ABSCOF (I, J)=TEMP 
IP(.NOT.OK) WRITE (6,530) I, J, TEMP 
IF (.HOT. OK) NEBROa=NERROR+1 

39 CONTINUE 

If DO is satisfied, thece were too many coeff. read. 

WRITE (6,540) KAEDON 
STOP 

All data cards have been read, check error count, W2ITS 
amount of data read. 

40 IP(HERfiOH.NE. 0) WRITE (6, 550) NEHROH 
IF (NEBROB.NE.O) STOP 
KABDS=KARDNO*KARDON-2 

WHITE (6,555) KABDS 

Form condensed abs. coeff. matrix. Allows all of data 
to be used. Each ACRED entry is the product of two . 
absorption coefficients. 

48 DO 70 J=1 .NGASES 

DO 60 K=J, NGASES 
SOW=0.0 

DO 50 1=1, LAMBDA 

SDM=SUM+ (ABSCOF (I,J) ♦ABSCOF (I, K) ) 


50 

CONTINUE 


ACRED (J,K) =SUM 


ACRED (K,J) =S0a 

60 

CONTINUE 
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70 


CONTINDE 


How find D(J) array. 

DO 90 J=1 ,NGASES 
SDI1=0.0 

DO 80 1=1, LAMBDA 

SaH=S0M> ABSCOF(I,J)*A5SOaP(I) 

80 CONTINUE 

D(J)=SUM 
90 CONTINUE 

Now trianqularize matrix to LL* form, 

K=0 

DO 130 J= UNGASES 
K=K*1 

TEMP=1/(DSQRT (ACHED (a,K))) 

KM1=K-1 

KDI7LM=NGASES-KM1 
DO 100 KDIV=1 ,KDIVLM 
KMDIV=KMUKDIV 

ACHED (J,KMDIV) =ACHED ( J, KMDIV) *TZMP 
100 CONTINUE 

D (J) =D (J) *TEMP 
LINDEX=J+1 
MINDEX=K<- 1 

DO 12C L=LINDEX, NGASES 
Q=-(ACRED(L,K) )/ACRED(J,K) 

DO 110 M=MIND£X, NGASES 

ACHED (L, M) = ACRED (L,M) +Q* (ACHED (J,M) ) 
no CONTINUE 

D(L)=D (L) >Q*D (J) 

120 CONTINUE 

130 CONTINUE 
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Zero entries in lower trianqle, all information is in upper 

DO 150 J=2,NGASES 
JH1=J-1 

DO 140 K=1,JM1 
ACHED (J, K) =0. 0 
140 CONTINOH 

150 CONTINOE 

c ■ ' 

C Now find solutions by back substitution, 

C Any neaative solutions are zeroea, and a record kept of 
C their existence. After all concentrations are found a 
C lump is made to a subroutone to remove all the absorp- 
C tion coefficients of these negative concentrations from 
C the elimination calculations. The exit is to a stage 
C where the ACRED array will be recalculated for a smaller 
C number of gases, 
c 

155 sua=o.o 
iacosc-0 

DO 170 J=1,NGAS£S 
JEI.Itt=NGASE3-J+1 

CONC (J£L,I.ai=iD^(>JELXHl nS!J.ai/ACE£.DX.J.£LJ.M 
'IP (CONC (JELIH) . LE.0.,0) IMCONC=IWCONC*1 
IFfCONC (JELIH) . LT.O. 0) CONC (JELIfl) =0 . 0 

sua=o.o 

JELial=JSLIM-1 
IFlJELiai.LI. 1) GO TO 17 1 
DO 160 N=JELia, NGASES 

SUa=SUM+ACEED(JELiai,N) ♦CCHC(N) 

160 CONTINUE 

170 CONTINUE 

171 XP1IMCONC) 172,172,200 
C 

C Now answers are in CONC array,. WRITE CONC in N/a*H and 
C torr, and give the ratios compared to the exp. estimates. 
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172 »BITE(6,560)NGASES 
DO 175 J=1,ilGASES 

WRITE(6,580)I8DEX(J) ,CONC(J) 

175 CONTINUE 

DO 180 J=1.NGAS£S 

CONC (J) =CONC(J) /133. 322 
180 CONTINUE 

HRITE(b, 570) NGASES 
DO 185 J=1,NGASES 

WEITE(t), 580) INDEX (J) , CONC (J) 

135 CONTINUE 
WRITE(6,590) 

1=1 

' DO 1 90 J=1 , NGASES 
READ (3,600) TEHP 
IP (INDEX (I) . NE.J) GO TO 190 
IF(TEaP. EQ.O. 0) GO TO 189 
TEMP=CCNC (J) /lEMP 
MSITE(6,610) INDEX (I) , TEMP 

189 1=1+1 ' 

190 CONTINUE 

STOP ^ 

This subroutine removes AESCOF of zeroed concentration 
from the calculations. Index is, the marKer showing 
which of the gases have been removed. 

200 NCONC=NGASSS 

201 DO 204 I=1,NCONC 

IP (CONC (I)) 205,205,204 

204 CONTINUE 
GO TO 48 

205 DO 210 J=I, NGASES 

JP1=J+1 

DO 208 K=l,LAflBDA - ~ 



AflSCDF (K, J) =4BSC0P (K, JP1) 

208 CONTINUE 

IMDSX (J) =INDBX (JP1) 

210 CONTINUE 

HGAS£S=NGASS3-1 
GO TO 204 

500 FOSilAT{2I2,2G6.3) 

510 PORHATCO* ,»Ea80R IN CONTROL PARAaETSES*) 

515 FOaaAT(2I2,G8.3) 

520 FORBAT i2I2,G10. 7) 

530 FORBAK IX, • ERROR IN CARE NITH 1= *,I2,' ,J= * , 12 , - 

** VALU3= » ,1PE14. 6) 

540 FORBATCO*, ‘TOO MUCH DATA INPUT’) 

550 FORBATCO', 'ERRORS FOUND IN ’,12,' DATA CARDS') 

555 FORBATC 1X,I3, ' DATA POINTS READ IN.') 

560 FORBATCO', ' THE CONCENTRATIONS OF THE ',12,' GASES 
♦ARE (IN N/B*M) '/) 

570 FORBATCO',' THE CONCENTRATIONS OF THE ',12,' GASES 
♦ARE(IN TORE ) ' /) 

580 P0RMAT(U,I2,5X,S14. 7) 

590 FORBAT( IX, 'THE RATIOS OP THE CALC, CQNC, 10 THE EXP 
♦2RIBSNTAL CONC. ARE;'/) 

600, FORBAT( 1X,G7.4) 

610 FORMAT (IX, 12, 5X, 314, 8) 

END 
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